arXiv:astro-ph/0512465v2 28 Jun 2006 


Draft version February 3, 2008 

Preprint typeset using D-T^;jX style emulateapj v. 6/22/04 


THE MASS ASSEMBLY HISTORY OF FIELD GALAXIES: DETECTION OF AN EVOLVING MASS LIMIT 

FOR STAR FORMING GALAXIES 

Kevin Bundy\ Richard S. Ellis\ Christopher J. Conselice\ James E. Taylor^ Michael C. Cooper^, 
Christopher N. A. Willmer^, Benjamin J. Weiner'^, Alison L. Coil^’®, Kai G. Noeske®, Peter R. M. Eisenhardt'" 

Draft version February 3, 2008 

ABSTRACT 

We characterize the mass-dependent evolution in a large sample of more than 8,000 galaxies using 
spectroscopic redshifts drawn from the DEEP2 Galaxy Redshift Survey in the range 0.4 < z < 1.4 and 
stellar masses calculated from AT-band photometry obtained at Palomar Observatory. This sample 
spans more than 1.5 square degrees in four independent fields. Using restframe {U — B) color and 
[OH] equivalent widths, we distinguish star-forming from passive populations in order to explore 
the nature of “downsizing”—a pattern in which the sites of active star formation shift from high 
mass galaxies at early times to lower mass systems at later epochs. Over the redshift range probed, 
we identify a mass limit, Mq, above which star formation appears to be quenched. The physical 
mechanisms responsible for downsizing can thus be empirically quantified by charting the evolution in 
this threshold mass. We find that Mq decreases with time by a factor of ~3 across the redshift range 
sampled according to Mq oc (1 -I- We demonstrate that this behavior is quite robust to the use 

of various indicators of star formation activity, including morphological type. To further constrain 
possible quenching mechanisms, we investigate how this downsizing signal depends on local galaxy 
environment using the projected 3’''^-nearest-neighbor statistic Dp 3 which is particularly well-suited 
for large spectroscopic samples. For the majority of galaxies in regions near the median density, there 
is no signihcant correlation between downsizing and environment. However, a trend is observed in the 
comparison between more extreme environments that are more than 3 times overdense or underdense 
relative to the median. Here, we find that downsizing is accelerated in overdense regions which host 
higher numbers of massive, early-type galaxies and fewer late-types as compared to the underdense 
regions. Our results signihcantly constrain recent suggestions for the origin of downsizing and indicate 
that the process for quenching star formation must, primarily, be internally driven. By quantifying 
both the time and density dependence of downsizing, our survey provides a valuable benchmark for 
galaxy models incorporating baryon physics. 

Subject headings: cosmology: observations, galaxies: formation, galaxies: evolution 


1. INTRODUCTION 

The redshift interval from z 1 to z = 0 accounts 
for roughly half of the age of the universe and pro¬ 
vides a valuable baseline over which to study the final 
stages of galaxy assembly. From many surveys spanning 
this redshift range, it is now well-established that the 
global star form ation rate (SFR) declines by an order of 
magnitude (e.g.,|B rQadIi!iiaL£LaLll^33lLillY_£LalJl22£ 

' de et allll 994lF!ore^tani^9^TwdsoT^tani2f ^ 

An interesting characteristic of this evolution in SFR is 
the fact that sites of active star formation shift from in¬ 
cluding high mass galaxies at z > 1 to only lower mass 
galaxies sub sequently. This patt ern of star formation, re¬ 
ferred to bv iGowie et a D iTT^ as “downsizing,” seems 


contrary to the precepts of hierarchical structure forma¬ 
tion and so understanding the physical processes that 
drive it is an important problem in galaxy formation. 

The observational evidence for the downsizing of star 
formation activity is now quite extensive. The pri¬ 
mary evidence comes from field surveys encompass- 
ing all classes of galaxies to z » 1 and beyond 


(IBrinchmann Rr. Ellid200fllBell et alJ2n051:i 

Bauer et al 

20051 l.lunea.u et aIJl2005l iFaber et aLll2005l 

Borch et al 


[a . Howeyer, the trends are also seen in specific popu¬ 


lations su ch as held spheroidals, both in their stell ar mass 
functions llFontana et alJ2004 lBundy etjilj2005|lajjdm 
more precise fundmnent al nlane studies ijTmT et al.l2005t 
lyan der Wei et al ] |20^ which track the eyolying mass- 
to-light ratio as a function of dynamical mass. The lat- 


Electronic address: kbundy@astro.caltech.edu, rse@atro.caltech.edu, 

^ Caltech MC 105—24, 1201 E. California Blvd., Pasadena, CA 
91125 

^ Department of Astronomy, University of California at Berke¬ 
ley, MC 3411, Berkeley, CA 94720 

® Steward Observatory, University of Arizona, Tucson, AZ 
98721 

^ Department of Astronomy, University of Maryland, College 
Park, MD 20742 

® Hubble Fellow, Steward Observatory, University of Arizona, 
Tucson, AZ 85721 

® UCO/Lick Observatory, University of California, Santa Cruz, 
CA 95064 

^ Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA 91109 


cc@feStiSlBSlt6SifH<flSl,IJ^®aM@(Spil6iifM.<(i&i^ CQSppiBh^tk'tSilCKblMifepIdu, cnaw@as. 

their star formation to within a few percent prior to z ~1, 
whereas lower mass ellipticals continue to grow by as 
much as 50% in terms of stellar mass after z ~ 1. Finally, 
detailed analyses of the spectra of nearby galaxies sug¬ 
gest similar trends llHeayens et al.ll2f)?)3: iJimenez et alJ 

l2f)nl . 

Downsizing is important to understand as it signifies 
the role that feedback plays in the mass-dependent eyolu- 
tion of galaxies. As a consequence, its physical origin has 
receiyed much attention theoreti cally. Recent analytic 
work by lDekel fc BirnboimI ll2004ll . for instance, suggests 
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that the distinction between star-forming and passive 
systems can be understood via several characteristic 
mass thresholds governed by the physics of clustering, 
shock heating and various feedback processes. Some 
of these processes have been implemented in numerical 
and semi-analytic mo dels, including ma ss-dependent 
star formation rates llMenci et al.l I200 !tI 1 . regulation 
throug^h feedback by s u pernovae (e.e.. IC ole et al.l 
120001 Benson et 1200, H INagashima fc YosM il 1200411 
and a ctive galactic nuclei ( AGN) (e.g.. ISilk fc Beef! 
1998LlGranatoetalJ_|200^_lDeke^_^_Birnboim|_|2004 

Hopkins et all 120051 ICroton et all 120051 IDe Lucia et all 

200!t 1 iScannapieco et al.l 1200,^ iBower et alJ 1200,4 

Cattaneo et all 1200011 . However, most models have, 
until now, primarily addressed the mass distinction 
between star-forming and qu iescent galaxies as define d 
at the present epoch (e.g., iKauffmann et alJ l200,Sbll . 
Quantitative observational measures of the evolving 
mass dependence via higher redshift data have not been 
available. 

This paper is concerned with undertaking a system¬ 
atic study of how the mass-dependent evolution of galax¬ 
ies progresses over a wide range of epochs. The goal 
is to quantify the patterns by which evolution proceeds 
as a basis for further modeling. Does downsizing result 
largely from the assembly history of massive early-type 
galaxies or is there a decline in the fraction of massive 
star-forming systems? In the quenching of star forma¬ 
tion, what are the primary processes responsible and how 
are they related to the hierarchical framework of struc¬ 
ture assembly as envisioned in the CDM paradigm? Does 
downsizing ultimately result from internal physical pro¬ 
cesses localized within galaxies such as star formation 
and AGN feedback, or is it caused by external effects 
related to the immediate environment, such as ram pres¬ 
sure stripping and encounters with nearby galaxies in 
groups and clusters? 

In this paper, we combine the large spectroscopic sam- 
ple contained in the DEEP2 Galaxy Redshift Survey 
ijPa.vis et filJl200^ with stellar masses based on exten¬ 
sive near-infrared imaging conducted at Palomar Obser¬ 
vatory to characterize the assembly history and evolution 
of galaxies since z ~ 1.2. Our primary goal is to quantify 
the downsizing signal in physical terms and test its envi¬ 
ronmental dependence so that it is possible to constrain 
the mechanisms responsible. A plan of the paper is as 
follows. Section m presents the observations and char¬ 
acterizes the sample while isi describes our methods for 
measuring stellar masses, star formation activity and en¬ 
vironmental density. We discuss how we estimate errors 
in the derived mass functions in 0 and present our re¬ 
sults in 0 We discuss our interpretations of the results 
in 0 and conclu de in 0 Whe re necessary, we assume 
a Chabrier IMF (|Chabrierll2flfl.‘lll and a standard cosmo¬ 
logical model with Om = 0.3, Ha = 0.7, Hq = lOOh 
km s“^ Mpc“^ and h — 0.7. 

2. OBSERVATIONS AND SAMPLE DESCRIPTION 

Constraining the processes that govern downsizing re¬ 
quires a precise measure of the evolving stellar mass func¬ 
tion of galaxies as a function of their physical state and 
environmental density. Achieving this ambitious goal re¬ 
quires multi-wavelength observations capable of reveal¬ 
ing quantities such as the ongoing star formation rate in 


a large enough cosmic volume to reliably probe a range 
of environments. Among these requirements, two ob¬ 
servational components are absolutely essential: a large 
spectroscopic survey and near-IR photometry. 

Spectroscopic redshifts not only precisely locate galax¬ 
ies in space and time, but enable the reliable determina¬ 
tion of restframe quantities such as color and luminosity. 
This in turn allows for accurate comparisons to stellar 
population templates which provide stellar mass-to-light 
ratios (M*/L) and the opportunity to conve rt from lumi¬ 
nosity to stellar mass. As demonstrated in iBundv et alJ 
(Unni), relying on photometric redshifts (photo-z) de¬ 
creases the typical precision of stellar mass estimates by 
more than a factor of three, with occasional catastrophic 
failures that lead to errors as large as an order of magni¬ 
tude. 

Spectroscopic redshifts are also c rucial for determinin g 
accurate environmental densities llCooner et al.l 1^0,41 . 
Even with the most optimistic photometric redshift un- 
certainties of = 0.02—COMBO-17 specifies « 0.03 
llWolf et alJ 1200,41 —a comparison between photo-z den¬ 
sity measurements and the real-space density in simu¬ 
lated data sets gives a Spearman ranked correlation co¬ 
efficient of only p = 0.4 (where p = 1 signifies perfect 
correlation, see iCooner et all 1200,4) . This uncertainty 
has the effect of smearing out the density signal in all 
but the lowest density environments in photo-z samples. 
For the spectroscopic DEEP2 sample used in this paper, 
Uz = 0.0001 and p = 0.8. 

Spectroscopic observations also provide line diagnos¬ 
tics that can discriminate the star formation activity 
occurring in galaxies. Given the various timescales in¬ 
volved, it is useful to identify ongoing or recent star 
formation activity by considering various independent 
methods including restframe (U—B) color, the equivalent 
width of [OH] A3727, and galaxy morphology. Gompar- 
isons between the different indicators highlight specific 
differences between early- and late-type populations de¬ 
fined in different ways. 

In addition to spectroscopy, the second essential ingre- 
die nt in this paper is near - IR ph otometry. As suggested 
byjKjuiffmajmj^_Ch ll 1 9981) and first exploited by 

iBrinchTnarmT^EIlifi (|2nnfij) . near-infrared and especially 
AT-band photometry traces the bulk of the established 
stellar populations and enables reliable stellar mass esti¬ 
mates for z < 1.5. With the addition of SED fitting from 
multi-band optical photometry, the uncertainty in such 
estimates can be reduced to factors of 2-3 based on Ri¬ 
band observations out to z « 1.5 llBrinchmann fc Elli^ 
[200^. The importance of AT-band observations is high¬ 
lighted for galaxies with z > 0.75, where the typical stel¬ 
lar mass uncertainty using the same technique applied to 
photometry with the z-band as the reddest filter would 
be worse by a factor of 3-4 (Bundy et ah, in prepara¬ 
tion). Uncertainties from optical mass estimates become 
even less secure as the redshift increases and increasingly 
bluer portions of the restframe SED are shifted into the 
reddest filter bands. Thus, the combined lack of AT-band 
photometry and spectroscopic redshifts leads to stellar 
mass errors greater than factors of 5-10 with catastrophic 
failures off by nearly two orders of magnitude. 

Because of these factors, the combination of the 
DEEP2 Galaxy Redshift Survey with panoramic IR 
imaging from Palomar Observatory represents the ideal 
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(and perhaps only) data set for tracing the evolution of 
accurately measured stellar masses and various indica¬ 
tors of star formation activity across different environ¬ 
ments to z ~ 1.4. We provide details on the specific 
components of this data set below. 

2.1. DEEP2 Spectroscopy and Photometry 

Th e DEEP 2 Galaxy Redshift Survey 
2003) utilizes the DEIMOS spectrograph 
20031) on the Keck-11 telescope and aims to measure 
~40,000 spectroscopic redshifts with z < 1.5 for galax¬ 
ies with Rab < 24.1. The survey samples four widely- 
separated regions, covering a total area of 3.5 square de¬ 
grees. Targeting of the spectroscopic sample was based 
on BRI photometry obtained at the Ganada-France- 
Hawaii Telescope (CFHT) w ith the 12Kx8K mosaic 
camera llCuillandre et al.ll200l|) . Gatalogs selected in the 
i?-band were constructed u sing the IMCAT photometric 
package llKaiser et al IM. and reac h a limiting magni¬ 
tude of j?.^R ^ 25.5 l|Goil et alJl2004D . The photometric 
calibration was computed with respect to Sloan Digital 
Sky Survey (SDSS) observations which overlap a portion 
of the GFHT sample. The observed colors, used to esti¬ 
mate the inferred restframe {U — B) colors in this paper, 
were measured using apertures defined by the object size 
in the i?-band image. Fu rther detai l s on t he GFHT pho¬ 
tometry can be found in iGoil et all 1)20041) . 

Sources were targeted for DEIMOS spectroscopy based 
on several criteria. As determined by magnitude, color, 
and size, objects in the photometric catalog were as¬ 
signed a probability, Pgai, for being a galaxy. Red- 
shift targets were required to have Pgai > 20% with 
magnitudes in the range, 18.5 < Rab < 24.1 (see 
IWillmer et alJl2r)0^ . In DEEP2 Fields 2-4, targets were 
also color-selected in {B — R) vs. {R — I) color space 
to lie predominantly at redshifts greater than 0.7. De¬ 
tails on the color cuts will be presented in Faber et al. 
(in preparation). These selection criteria successfully re¬ 
cover 97% of the Rab < 24.1 population at z > 0.75 with 
only "^10% contam ination from lower redshift galaxies 
ijPavis et alJl200^ . The redshift survey in these three 
fields is now complete, providing a total of 21592 suc¬ 
cessful redshifts over 3 square degrees. The last field, 
the Extended Groth Strip (EGS) covers 0.5 square de¬ 
grees and is currently 75% complete with a total of 9501 
galaxies in the range 0 < z < 1.5. The selection of targets 
in EGS is magnitude-limited, providing a valuable sam¬ 
ple for verifying the success of the color selection used in 
Fields 2-4. The sampling rate of galaxies satisfying the 
target criteria is ~60% and DEEP2 galaxies from all 4 
fields are included in this paper. 

DEEP2 redshifts were determined by comparing vari¬ 
ous templates to observed spectra as well as fitting spe¬ 
cific spectral features. This process is interactive and will 
be described in Newman et al. (in preparation). Spectro¬ 
scopic redshifts are used in this paper only when two or 
more features have been identi fied in a given sp ectrum 
(giving “zquality” values > 3. iFaber et alJl2?)fi^ . The 
fraction of objects for which this process fails to give a 
reliable redshift is roughly 30% and is dominated by faint 
blue galaxies, the majority of which are beyond z ~ 1.4 
where the [OH] A3727 feature is r edshifted beyond th e 
DEEP2 spectral wavelength range llWillmer et alJl^005D . 
More details on the redshift success rate are provided in 


IWillmer et alJ ll2005D . and IGoil et al.l ll2004D further dis¬ 
cuss the survey strategy and spectroscopic observations. 

2.2. Palomar Near-IR Imaging 

Motivated largely by this study, we have conducted an 
extensive panoramic imaging survey of all four DEEP2 
fields with the W ide Field Infrared Camera (WIRC, 
iWilson et al.ll200^ on the 5m Hale Telescope at Palo¬ 
mar Observatory. We describe the salient points of this 
survey here and will provide further observational details 
in a later article (Bundy et al., in preparation). 

The Palomar survey commenced in fall 2002 and was 
completed after 65 nights of observing in October 2005. 
Using contiguously spaced pointings (each with a cam¬ 
era field-of-view of 8f6x8f6) tiled in a 3x5 pattern, we 
mapped the central third of Fields 2-4 to a median 80% 
completeness depth greater than Kab = 21.5, with 5 
pointings deeper than Kab = 22.5. The imaging in 
Fields 2-4 accounts for 0.9 square degrees or 55% of the 
Palomar As-band survey. 

The remainder of the data was taken in the EGS where 
the Ag-band data covers 0.7 square degrees, but to var¬ 
ious depths. The EGS was considered the highest pri¬ 
ority field in view of the many ancillary observations— 
including HST, Spitzer, and X-ray imaging—obtained 
there. As the orientation of the WIRC camera is fixed on 
the sky and the EGS traces a ~45 degree strip 16' wide, 
to fully cover the spectroscopic field in the E-W direc¬ 
tion requires rows of three WIRC pointings. The N-S 
direction requires about 12 different positions, so, in to¬ 
tal, 35 WIRC pointings were used to map the EGS in 
the Ag-band. The deepest observations were obtained 
along the center of the strip where there is complete over¬ 
lap between WIRC and the spectroscopically-surveyed 
area. In these regions, the typical depth is greater than 
Kab = 22.5. The rest of the southern half of the EGS 
reaches Kab = 22.3, while that for the northern half is 
complete to Kab = 21.7. Additional Palomar J-band 
observations were obtained for most of the central strip 
of the EGS and for Fields 3 and 4. These provide J-band 
photometry for roughly half of the Ag-band sample and 
are useful in improving photometric redshifts. 

Individual Ag-band exposures at a given pointing were 
taken with an integration time of 2 minutes using coaddi¬ 
tions of either 4 x 30 seconds in average temperature con¬ 
ditions, 3 frames x 40 seconds in cooler conditions, or 6 
frames x 20 seconds in warm conditions. The exposures 
were dithered at 27 positions in a non-repeating, pseudo¬ 
random «7"pattern and combined into 54-minute mo¬ 
saics using a double-pass reduction pipeline developed 
specifically for WIRC. At a given pointing, individual 
mosaics were often obtained on different nights and so 
may vary in terms of seeing, sky background levels and 
transparency. Most Ag-band pointings consist of more 
than two independently-combined mosaics with the deep¬ 
est pointings comprising as many as 6 independent mo¬ 
saics. Mosaic coaddition was performed using an algo¬ 
rithm that optimizes the depth of the final image by 
applying weights based on the seeing, background, and 
transparency of the subframes. The final seeing FWHM 
in the Ag-band images ranges from 0'.'8 to 1'.'2, and is 
typically I'.'O. Photometric calibration was carried out 
by referencing standard stars during photometric condi¬ 
tions and astrometric registration was performed with 


Davis_££^ 

Faber et al. 
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respect to DEEP2 astrometry Csee lCoil et al.ll20ti^ us¬ 
ing bright stars from the CEHT i?-selected catalog. More 
details on the survey strategy and image processing are 
presented in Bundy et al. (in preparation). 

After masking out the low signal-to-noise perime- 
ter of the final A«-ban d images, we used SExtractor 
llBertin fc Arnoiit^ll 99611 to detect and measure the K^- 
band sources. As an estimate of the total Ag-band mag¬ 
nitudes, which are used to derive the luminosity and 
stellar mass of individual galaxies (see we use the 

MAG_AUTO output from SExtractor. We do not ad¬ 
just this Kron-like magnitude to account for missing light 
in extended sources, which introduces two slight biases 
in our total magnitudes. Eirst, the average difference 
between the MAG_AUTO values of a given galaxy and 
its corresponding magnitude measured in a 4" diameter 
aperture is wO.OS mag fainter at the highest redshifts 
compared to the lowest redshifts in the sample. This 
effect leads to a systematic underestimate of 0.01 dex 
in stellar masses at high -2 compared to low-z. In ad¬ 
dition, MAG_AUTO systematically underestimates the 
total magnitude of b lue galaxies (defined by restframe 
U — B color, see in the lowest redshift interval as 
compared to red galaxies. The magnitude of this effect, 
which likely results from the more extended nature of 
typical blue galaxies, is also wO.OS mag but is not signifi¬ 
cant in the middle or high redshift intervals. Considering 
these systematic effects contribute less than 0.01 dex to 
the final stellar mass estimates, we do not explicitly cor¬ 
rect for them. 

In addition to total magnitudes, we also use SExtrac¬ 
tor to measure aperture photometry in diameters of 2", 
3", 4", and 5". To determine the corresponding magni¬ 
tudes of A^-band sources in the CEHT BRI and Palomar 
J-band images, we applied the IDL photometry proce¬ 
dure, APER, to these data, placing apertures with the 
same set of diameters at sky positions determined by 
the Ag-band detections (all magnitudes were corrected 
to airmass = 1). This method was adopted with the 
aim of producing a A-selected catalog—about 25% of 
the Ag-band sources do not have optical counterparts in 
the CEHT IMCAT catalog. Magnitudes measured with 
APER are consistent with SExtractor aperture magni¬ 
tudes, and experimentation with color-color diagrams 
and photometric redshifts demonstrated that the 2" di¬ 
ameter colors exhibited the least scatter. We therefore 
use the 2 " aperture colors for fitting template spectral 
energy distributions (SEDs) to constrain stellar mass es¬ 
timates, as well as to estimate photo-z’s for sources be¬ 
yond the DEEP2 spectroscopic limit. 

Photometric errors and the Ag-band detection limit of 
each image were estimated by randomly inserting fake 
sources of known magnitude into each Ag-band image 
and recovering them with the same detection parameters 
used for real objects. The inserted objects were given 
Gaussian profiles with a EWHM of I'.'S to approximate 
the shape of slightly extended, distant galaxies. We de¬ 
fine the detection limit as the magnitude at which more 
than 80% of the simulated sources are detected by SEx¬ 
tractor. Robust photometric errors based on simulations 
involving thousands of fake sources were also determined 
for the BRI and J-band data. These errors are used to 
determine the uncertainty of the stellar masses and in the 
determination of photometric redshifts, where required. 


2.3. The Primary Sample 

Given the various ingredients necessary for the data 
set, it is helpful to construct separate samples based on 
the differing completeness limits for the A-limited spec¬ 
troscopic sample and the A-limited Palomar catalog. 

We will define the primary sample as that compris¬ 
ing galaxies with DEEP2 spectroscopic redshifts that are 
also detected in the Palomar Ag-band imaging. Red¬ 
shift counterparts were found by cross-referencing the 
Ag-band catalog with the DEEP2 redshift catalog. A 
conservative tolerance of 1"1 in the difference between 
the Ag-band and DEEP2 positions was used to select 
matches, although the offset for most sources was less 
than 0'.'5. The relatively low surface density of both cat¬ 
alogs assures that the number of spurious associations 
is less than 1-2 percent. The fraction of DEEP2 red¬ 
shifts detected in the Ag-band varies from «65% for Ag- 
band depths near Kab = 21.7 to R:!90% for Kab = 22.7. 
After removing the redshift survey boundaries to allow 
for unbiased density estimates (see ^3.311 . Fields 2, 3, 
and 4 host 953, 1168, and 1704 sample galaxies respec¬ 
tively, all with secure spectroscopic redshifts in the range 
0.75 < 2 < 1.4. The fourth field, the EGS, contains 4770 
galaxies with redshifts in the range 0.2 < 2 < 1.4. 

In the analysis that follows, we divide this sample into 
three broad redshift intervals. The first, 0.4 < 2 < 0.7, 
contains 943 galaxies, drawn entirely from the EGS while 
the second (0.75 < 2 < 1.0, 2210 galaxies) and third 
(1.0 < 2 < 1.4, 1430 galaxies) draw from all four fields 
(see Table P). Although the entire spectroscopic sam¬ 
ple was selected to have Rab < 24.1, the limiting Ag- 
band magnitude was chosen separately for each redshift 
interval. Because the Palomar Ag-band survey covers 
different areas to different depths, the area and vol¬ 
ume of a given subsample depend on the Ag-band limit 
that is applied. This is advantageous and allows us to 
choose limits for each redshift bin that maintain adequate 
statistics and stellar mass completeness while producing 
subsamples that probe similar cosmological volumes—an 
important consideration for making environmental com¬ 
parisons in a consistent fashion. In the three redshift 
bins, we select galaxies with secured Ag-band detections 
brighter than 21.8, 22.0, and 22.2 (AB), respectively. In 
the standard cosmology we have adopted with h — 0.7, 
the areas sampled with these limits result in comoving 
volumes of 0.5, 1.4, and 2.3 in units of 10® Mpc®. 

2.4. The Photo-z Supplemented Sample 

As mentioned previously, photometric redshifts 
are insufficient fo r accurate density measurements 
l|CooDer et al.ll200^ . do not offer [OII]-based SFR esti¬ 
mates j_and_signi&ant]y decrease the precision of stellar 
mass llBundv et iOlMl) and restframe color estimates. 
However, photometric redshifts do offer the opportunity 
to augment the primary sample and extend it to fainter 
magnitudes, providing a way to test for the effects of 
incompleteness in the primary spectroscopic sample be¬ 
cause of the various magnitude limits and selection pro¬ 
cedures. With this goal in mind, we constructed a com¬ 
parison sample using the optical-l-near-IR photometry 
described above to estimate both photometric redshifts 
and stellar masses. We discuss the importance of this 
comparison in interpreting our results in section 
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TABLE 1 
Sample statistics 



Primary Spec- 2 :, Rab ^ 24.1 

Photo- 2 : Supplemented, Rab ^ 25.1 

Sample 

Number 

Number 

fspec-z 

/aNNz 

/bpz 

EGS Field 

Full sample 

2669 

8540 

36% 

51% 

13% 

0.4 < 2 < 0.7 

943 

3026 

36% 

62% 

2% 

0.75 < 2 < 1.0 

1003 

2801 

42% 

46% 

12% 

1.0 < 2 < 1.4 

723 

2713 

29% 

43% 

28% 


DEEP2 Fields 2 

-4 




Full sample 

1914 

10156 

21% 

68% 

11% 

0.4 < 2 < 0.7 

— 

4264 

2% 

97% 

1% 

0.75 < 2 < 1.0 

1207 

2865 

42% 

48% 

10% 

1.0 < 2 < 1.4 

707 

3027 

27% 

47% 

26% 


Note. — The listed values reflect cropping th e DE EP2 survey borders to allow for ac¬ 
curate environmental density measurements (see and, for the three redshift intervals, 
As-band magnitude cuts of 21.8, 22.0, and 22.2 (AB). 


Photometric redshifts were derived in two ways. First, 
because the DEEP2 multi-slit masks do not target ev¬ 
ery available galaxy, there is a substantial number of ob¬ 
jects that satisfy the photometric criterion of Rab < 24.1 
without spectroscopic redshifts. These are ideal for neu¬ 
ral network photo-z estimates because of the availability 
of a large spectros copic training set. For these galax¬ 
ies, we use ANNz llCollister fc Laha'ill2004ll to measure 
photo-z’s, training the network with the EGS spectro¬ 
scopic sample so that the ANNz results cover the full 
range, 0.2 < z < 1.4, without being subject to biases 
from any color selection. Based on comparisons to the 
spectroscopic samples in Fields 2-4, the ANNz results 
are in excellent agreement with Az/(1 -|- z) « 0.07. 

While the Rab < 24.1 ANNz photo-z results present a 
complete sample, they do not contain fainter galaxies be¬ 
cause no spectroscopic sample is available to train them. 
We therefore define a second sample of galaxies with 
24.1 < Rab < 25.1 and 3 (t detections in th e BIK bands . 
For these we use the BPZ photo-z package llBeni tezl2000ll 
and 2'.'0 diameter aperture photometry, including J-band 
where available. We first optimize BPZ by comparing to 
spectroscopic samples. These tests reveal a systematic 
offset in the BPZ results that likely ar ises from the as¬ 
sumed HDF-North prior llBenftedl2000ll . We fit a linear 
relation to this offset and remove it from all subsequent 
BPZ estimates. This improves the spectroscopic compar¬ 
ison to Az/(l-|-z) « 0.17. For sources with Rab < 24.1, 
we find good agreement between the ANNz and BPZ re¬ 
sults, although the BPZ estimates exhibit a larger degree 
of scatter. 

Using the ANNz and BPZ results, we construct photo- 
z supplemented samples in each redshift bin with an R- 
band limit of Rab < 25.1. The first bin, with 0.4 < 
z < 0.7, contains 7290 galaxies; 16% have spectroscopic 
redshifts (available in the EGS only), 83% have ANNz 
photo-z’s, and only 1% have BPZ photo-z’s. The second 
bin, with 0.75 < z < 1.0, contains 5666 galaxies; 42% 
have spectroscopic redshifts, 48% have ANNz photo-z’s, 
and 10% have BPZ photo-z’s. The third bin, with 0.75 < 
z < 1.0, contains 5740 galaxies; 28% have spectroscopic 


redshifts, 45% have ANNz photo-z’s, and 27% have BPZ 
photo-z’s. The sample statistics are summarized in Table 

m 

3. DETERMINING PHYSICAL PROPERTIES 

The goal of this paper is to derive key physical quan¬ 
tities that can be used to describe the stellar mass, evo¬ 
lutionary state, and environmental density of galaxies, 
and to use these measures to understand the physical 
processes that drive the broad patterns observed. In this 
section we discuss our methods for reliably determining 
these key variables by making use of the unique combi¬ 
nation of spectroscopy and near-IR imaging offered by 
the DEEP2/Palomar survey described above. The un¬ 
certainties discussed below refer to the primary sample. 
Our full error analysis is described in SH 

3.1. Stellar Masses Estimates 


AT-band luminosities alone provide stell ar mass esti¬ 
mate s that are uncertain by factors of 5-7 llBrinchmannI 
[1993). However, for samples of known spectroscopic 
redshift, using optical-infrared color information can 
constrain the stellar population and M^/L ratio so as 
to reduce this uncertainty to ~0.2 dex. In this pa¬ 
per, we e mploy such a t echni que using the method de¬ 
scribed in iBundv et all ll2005li and discussed further in 
Bundy et al. (in preparation). We use a code that is 


Ba yesian in nature a n d based on th e precepts described 
in jl^iiffmann et abl (|2nn3all and iBrinchmann fc Elliij 

Briefly, the code uses BRIK colors (measured using 
2'.'0 diameter aperture photometry matched to the K- 
band detections) and spectroscopic redshifts to compare 
the observed SED of a sample galaxy to a grid of 13440 
synthetic SEDs Ifrom lBriizual fc Gharlotl2003D spanning 
a range of star formation histories (parametrized as an 
exponential), ages (restricted to be less than the age of 
the uniyerse at the obseryed redshift), metallicities, and 
dust content. The choice of models and population syn¬ 
thesis code emplo yed can introduce system atic uncertain¬ 
ties. For example. iMarastonet 'n]0ool show that p op- 
ular recipes such as those in lBruziial fc Gharloil i2003H do 







































6 


Bundy et al. 


not accurately chart the evolution of Thermally-Pulsing 
Asymptotic Giant Branch stars. For young template 
models (with ages less than 2 Gyr), the Maraston models 
predict stellar masses that can be lower by as much as 
^60%. We note this potential source of error, but do not 
explicitly correct for it in estimating stellar masses here. 

Once the grid of models has been defined, the Kg- 
band M^^/Lk, minimum and the probability that 
each model accurately describes a given galaxy is calcu¬ 
lated at each grid point. The corresponding stellar mass 
is then determined by scaling ratios to the Kg- 

band luminosity based on the total iO-band magnitude 
and redshift of the observed galaxy. The probabilities are 
then summed (marginalized) across the grid and binned 
by model stellar mass, yielding a stellar mass probability 
distribution for each sample galaxy. We use the median 
of the distribution as the best estimate. 

The stellar mass measured in this way is robust to de¬ 
generacies in the model, such as those between age and 
metallicity. Although these degeneracies can produce bi- 
modal probability distributions, even in these cases, the 
typical width of the distribution gives uncertainties from 
the model fitting alone of 0.1-0.2 dex. For about 2-3% 
of the SED fits, the minimum values are significantly 
greater than 1.0. For these more poorly-constrained ob¬ 
jects, we add an additional 0.2 dex in quadrature to the 
final mass uncertainty. Although in principle, the best 
fitting model also provides estimates of the age, metallic¬ 
ity, star formation history, and dust content of a sample 
galaxy, these quantities are much more affected by de¬ 
generacies and are poorly constrained compared to the 
stellar mass. 

We also note that real galaxies are likely to have more 
complex star formation histories than the simple stellar 
populations assumed here. However, we do not attempt 
to fit the observed SEDs with more complicated models, 
including those with multiple components and bursts, be¬ 
cause our photometric observations do not provide suffi¬ 
cient constraints. More complex models are particularly 
relevant to galaxies with recent bursts of star formation 
which can hide underlying stellar mass. Several authors 
have investi gated this effect. Stud ying the most extreme 
possibilities IShaolev et IDiEonS), for example, find that 
the stellar mass of blue galaxies with recent starbursts 
could be underestimated by a fact or of 5. More typ ical 
values are probably much less, and iTran et al.l ll2004ll es¬ 
timate the fraction of post starbursts identified as E-l-A 
galaxies in t he range 0.3 < z < 1 to be a few percent 
(although see lYan et alJl200,^ . so the effect of starbursts 
in our sample is likely to be small. 

Photometric errors enter the stellar mass analysis by 
determining how well the template SEDs can be con¬ 
strained by the data. Larger photometric uncertainties 
“smear out” the portion of the model grid that fits an ob¬ 
served galaxy with high probability. This is reflected in 
a wider stellar mass probability distribution. Additional 
uncertainties in the iG-band luminosity (from errors in 
the observed total Ag-band magnitude) lead to final stel¬ 
lar mass estimates that are typically good to 0.2-0.3 dex. 
The largest systematic source of error comes from the 
assuin ed IMF, in this case that proposed by iGhabrieJ 
iMoS). The stellar masses we derive using this IMF can 
be converted to Salpeter by adding ~0.25 dex, although 
we note this offset represents an average with a scatter 



(U-B),.,. 



Fig. 1. — The restframe {U — B) color distribution in the middle 
and high redshift bins. Overplotted are the distributions of galax¬ 
ies whose star formation rate as derived from the [Oil] equivalent 
width is greater than 0.2 MQ/yr (blue dashed line) and less than 
0.2 M 0 /yr (red dotted line). The (U — B) bimodality discriminant 
and its Itr range is indicated by the vertical dashed line and grey 
shading. 

across the models in our grid of ~0.06 dex. 

3.2. Indicators of Star Formation Activity 

In this paper, we adopt two independent approaches 
for identifying those galaxies that are undergoing, or have 
recently experienced, active star formation. The first is 
the restframe (U — B) color, estima ted with the same 
methods as in iWillmer et al.l 1)200,^ and based on the 
IMCAT photometry measured for the GFHT BRI data. 
The fc-corrections, which translate observed BRI colors 
into restframe (U—B) values, are determin ed by compar¬ 
ison to a set of 43 nearby galaxy SEDs from IKinnev et aP 
lIT^ . Second-order polynomials were used to estimate 
(U—B) and fc-c orrections as a f unctio n of observed color. 
Appendix A in IWillmer et al.l l|200,^ provides more de¬ 
tails on this technique. 

The {U — B) color distribution exhibits a clear bi¬ 
modality that is used to divide our sample into “Blue” 
(late-type) and “Red” (early-type) subsamples using 
the same luminosi t y-dep endent color cut employed by 
Ivan Dokkiim et H] feoOl^ : 

U-B = -0.032(Mb -b 21.52) -b 0.454 - 0.25 (1) 

This formula, defined in the Vega magnitude system, di¬ 
vides the sample well at all redshifts (see Figure QJ, and 
so we do not apply a co rrection for potential rest-frame 
evolution with redshift (IWillmer et al.ll2005D . We note 
that the proportion of red galaxies at high-z in Figure 
n]is in part suppressed by the i?-band selection limit of 
Rab = 24.1 (selection effects are discussed in detail in 

For galaxies with z > 0.75, the [Oil] 3727 A emission 
line falls within the spectroscopic range probed by the 
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DEEP2 Survey and this can be used to provide a sec¬ 
ond, independent estimate of the SFR. For these galax¬ 
ies, we measure the intensity of the [Oil] emission line 
by fitting a double gaussian—with the wavelength ratio 
constrained to that of the [Oil] doublet—using a non¬ 
linear least squares minimization. We measure a robust 
estimate of the continuum by taking the biweight of the 
spectra in two wavelength windows each 80 A long and 
separated from the emission line by 20 A. We consider 
only spectra where the equivalent width of the feature is 
detected with greater than 3cr confidence. 

Because the DEEP2 spectra are not flux calibrated, we 
use the following formula from iGuzmau et a, 1.1 ill OOTTi to 
estimate the [Oil] SFR: 

SFR(Mq yr-i) « 1O-“'6-° '‘(*^^-^«0 )EW[oii] (2) 

This relation utilizes the restframe Mn estimated in the 
same way as the {U — B) colors I'see lWillmer et al]l200fi|) 
and provides an approximate value for the SFR with¬ 
out correcting for metallicity effects, which can intro¬ 
duce random and systematic error s of at least 0.3-0. 5 
dex in the SFR derived from L[oii] llKewlev et al.ll200'^ . 
In addition, Equation [21 was opti mized for a sample o f 
compact emission line galaxies Isee lGuzman et id]ll997l) 
which likely differs in the amount of extinction and typ¬ 
ical [Oil]/Ha ratio compared to the sample here, yield¬ 
ing SFR estimates that could be SYStem atically off by a 
factor of ^3. Moreoyer, recent work by IPanovich et ahl 
ll20f)l demonstrates that SFR’s based on re-radiated IR 
radiation can be orders of magnitude larger than opti- 
cal/UV estimates, and studies of the AGN population 
in DEEP2 indicate that [OH] emission is often associ¬ 
ated w ith AGN, furth er biasing [OH]-based SFR esti¬ 
mates llYan et al.ll2005ll . Future efforts, especially those 
utilizing the multi-wayelength data set in the EGS, will 
be useful in refining the SFR estimates. In the present 
paper, we accept the limitations of Equation (21 for de¬ 
termining absolute SFR’s, noting that our primary re¬ 
quirement is not a precision estimate of the SFR for each 
galaxy but only the broad diyision of the field population 
into active and quiescent subsets. 

In support of this last point, Figure ^ compares the 
{U — B) and [OH] star formation indicators directly in 
the two higher redshift bins where both diagnostics are 
available. The solid histogram traces the full {U — B) 
distribution with the vertical dashed line (and shading) 
indicating the median and Icr range of the color bimodal¬ 
ity discriminant (the range results from the dependence 
on Mb in Equation^ . Using the independent diagnostic 
of the [OH] equivalent width, we can similarly divide the 
population into high (blue dashed histogram) and low 
(red dotted histogram) star-forming populations using a 
cut of SFR[oii] = 0.2 Mq yr“^, which is the median SFR 
of galaxies with 0.75 < z < 1.0. Figure Q] clearly shows 
the effectiveness of the {U — B) color cut in distinguish¬ 
ing the populations in both cases. The fraction of red 
galaxies in the high SFRpuj population is less than 8% 
in the middle redshift bin and less than 17% in the high- 
z bin. The fact that a non-zero fraction of blue galaxies 
is made up of the low SFRjqii] population is likely in¬ 
dicative of the 1-2 Gyr timescale required for galaxies 
to join the red sequence and suggests that even minor 
episodes of star formation can lead to blue restframe col- 



Log (S 3 / <E 3 >) 

Fig. 2.— The distribution in the relative overdensity in log units 
as measured across the sample by the 3^ ^-nearest-neighbor statistic 
as implemented bv ICooper et all 1200511 . The vertical dotted lines 
at ±0.5 dex divide the distribution into low, middle, and high 
density regimes (see text for details). 

ors llGebhardt et alJ2003r . In addition, the median value 
of the measured [OH] SFR increases in the high-z bin. 
This means that galaxies satisfying the SFR[oii] = 0.2 
Mq yr“^ cut at high redshift are more vigorously form¬ 
ing stars, and therefore bluer. 

Further details on the evolving SFR will be presented 
in Noeske et al. (in preparation). We also note that the 
iU—B) and [OH] star formation indicators are consistent 
with the star formation histories recovered by the SED 
template fitting procedure used to refine stellar mass es¬ 
timates and described in MJi This agreement is ex¬ 
pected because the restframe color and SED fitting are 
both determined by the observed colors. Such consis¬ 
tency demonstrates that the measured physical proper¬ 
ties that we use to divide the galaxy population are also 
reflected in the best-fit SED templates that determine 
stellar mass (see Q. 

3.3. Environmental Density 

Charting galaxy evolution over a range of environments 
represents a key step forward that can only be achieved 
through large spectr oscopic redshi f t surv eys such as in 
the DEEP2 survey. iCooper et alJ ll2005ll rigorously in¬ 
vestigate the question of how to provide precise environ¬ 
mental density estimates in the context of redshift sur¬ 
veys at z ~ 1. That work clearly shows, via comparisons 
to simulated samples, that large samples with photo-z’s 
are very poorly suited to providing accurate density mea¬ 
sures in all bu t the most underdens e environments. 

Specifically, iCooper et al.1 ll2005ll find that for envi- 
ronmental studies at high-z, the projected n*^-nearest- 
neighbor distance, offers the highest accuracy over 

the greatest range of environments. This measure is the 
field counterpart to projected density estimates first ap¬ 
plied to studies of cluster environments (e.g., iDressled 
HMg). The statistic is defined within a velocity window 
(Az) = ±1000 km s“^) used to exclude contaminating 
foreground and background galaxies. It is therefore par¬ 
ticularly robust to redshift space distortions in high den¬ 
sity environments without suffering in accuracy in un¬ 
derdense environments. In addition, the effect of survey 
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boundaries on Dp^n is easily understood and mitigated 
by excluding a small strip around survey edges. 

In this paper, we utilize the projected 3'’'^-nearest- 
neighbor distance, Dp^^^ excluding galaxies closer than 
1 h~^ Mpc (~3') from a survey boundary. The choice 
of n = 3 does not significantly affect Dp n which has a 
weak dependence on n in both high and low density en¬ 
vironments for n < 5 llCooDer et ai.ll200'^ . The effects of 
survey target selection must be carefully considered be- 
cause they can intro duce biases as a function of redshift. 
iCooper et al.1 (1200511 find that the sampling rate and se¬ 
lection function in the DEEP2 survey equally probes all 
environments at z ~ 1 in a uniform fashion. Although 
the DEEP2 survey secures redshifts for roughly 50% of 
galaxies at 2 ; ~ 1 , its sparse selection algorithm does 
not introduce a significant environment-dependent bias. 
However, the absolute value of Dp^^ will increasingly un¬ 
derestimate the true density with increasing redshift as 
the sampling rate declines. To handle this effect, we first 
convert Dp^^ into a projected surface density, E 3 , using 
^3 = ^/(TriDp g). We then calculate the relative over¬ 
density at the location of each galaxy by dividing the 
observed value of E 3 by the median surface density cal¬ 
culated in bins of Az = 0.04. The relative overdensity is 
thus insensitive to the DEEP2 sampling rate, providing 
a reliable statistic that can be compared across th e full 
redshift range of the sample l|Cooper et al.ll2fiC)^ . We 
note that a small bias in this statistic may be introduced 
by the fact that DEEP2 preferentially selects blue galax¬ 
ies as the redshift increases. Tests indicate this effect is 
smaller than the typical density uncertainty (a factor of 
^3) and is not likely to be important given the coarse 
density bins we use to divide the sam ple. Further details 
are provided in lCooner et alJ l| 200 (i(l . 

The distribution of the relative overdensity for our pri¬ 
mary sample is plotted in Figure [21 In the analysis that 
follows, we consider two ways of dividing the sample by 
environmental density. In the first, we separate galaxies 
according to whether they lie in regions above or below 
the median density (corresponding to a measured over¬ 
density of zero in Figure!^. This is the simplest criterion 
but does mean the bulk of the signal is coming from re¬ 
gions which are not too dissimilar in their environs. In 
the second approach, we divide the density sample into 
three bins as indicated by the vertical dotted lines in 
Figure (21 The thresholds of 0.5 dex, or 0.77cr, above 
and below the median density were chosen to select the 
extreme ends of the density distribution where the field 
sample begins to probe group and void-like environments 
that are 3-100 times more or less dense than average. 
The conclusions presented in are not sensitive to the 
precise location of these three thresholds. 

4. CONSTRUCTING THE GALAXY STELLAR MASS 
FUNCTION 

The DEEP2/Palomar survey presents a unique data 
set for constraining the galaxy stellar mass function at 
z ~ 1. Previous efforts have so far relied on smaller and 
more limited samples, often without spectroscopic red- 
shifts. Only via the combination of size, spectroscopic 
completeness, and depth, does the current sample en¬ 
able us to study the evolving relationships between stel¬ 
lar mass, star formation activity, and local environment 
in a statistically robust way. Our primary tool in this ef¬ 


fort is charting the galaxy stellar mass function of various 
populations. 

Deriving the stellar mass function in a magnitude lim¬ 
ited survey requires corrections for the fact that faint 
galaxies are not detected thr oughout the en tire survey 
volume. The Vmax formalism l|Schmidtll1 is the sim¬ 
plest technique for handling this problem but does not 
account for density inhomogeneities that can bias the 
shape of the derived mass function. While ot her tech¬ 
nique s address this problem (for a review, see iWillmeJ 
[1993), they face different challenges such as determining 
the total normalization. For sufficiently large samples 
over significant cosmological volumes, density inhomo¬ 
geneities cancel out and the Vmax method produces reli¬ 
able results. Considering this and the fact that we wish 
to use our data set to explicitly test for the effects of 
density inhomogeneities, we adopt the Vmax approach, 
which offers a simple way to account for both the R and 
iC-band limits of our sample. For each galaxy i in the 
redshift interval j, the value of V^ax is given by the min¬ 
imum redshift at which the galaxy would leave the sam¬ 
ple, becoming too faint for either the R or Ag-band limit. 
Formally, we define 

f^high 

= dn, —dz (3) 

where dVlj is the solid angle subtended by the sample 

defined by the limiting Ag-band magnitude, for 

the redshift interval j, and dV/dz is the comoving volume 
element. The redshift limits of the integral are: 

Zhigh = (4) 

Ziow = zi,,^ (5) 

where the redshift interval, j, is defined by [Zmin^ ^max]^ 
^’sfers to the redshift at which the galaxy would 
still be detected below the Ag-band limit for that partic¬ 
ular redshift interval, and zr,^^ is the redshift at which 
the galaxy would no longer satisfy the A-band limit of 
Rab < 24.1. We use the SED template fits found by 
the stellar mass estimator to calculate and ZR^.^, 

thereby accounting for the fc-corrections necessary to 
compute accurate Vmax values (no evolutionary correc¬ 
tion is applied). 

In constructing the Vmax mass functions, we also 
weight the spectroscopic sample to account for incom¬ 
pleteness resulting from the DEEP2 color selection and 
redshift s uccess rate. We close ly follow the method de¬ 
scribed in IWillmer et alJ l|2005ll . but add an extra dimen¬ 
sion to the reference data cube which stores the num¬ 
ber of objects with a given Ag-band magnitude. Thus, 
for each galaxy i we count the number of objects from 
the photometric catalog sharing the same bin in the 
{B — R)/{R — /) /Rab /Ag-band parameter space as well 
as the fraction of these with attempted and successful 
redshifts. As mentioned in EH «30% of attempted 
DEEP2 redshifts are unsuccessful, mostly because of 
faint, blue galaxies bey ond the redshift limi t accessible to 
DEEP2 spectroscopy ijWillmer et alJl2fln5fl . Failed red¬ 
shifts for red galaxies are more likely to be within the ac¬ 
cessible redshift range but simply lacking in strong, iden¬ 
tifiable features. We therefore use the “optimal” weight- 
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ing model llWillmer et al.ll200!TD . which accounts for the 
redshift success rate by assuming that failed redshifts 
of red galaxies (defined by the {U — B) color bimodal¬ 
ity) follow the same distribution as successful ones, while 
blue galaxies with failed redshifts lie beyond the redshift 
limit {z « 1.5) of the sample. The final weights are then 
based on the probability that a successful redshift would 
be obtained for a given galaxy. They also account for the 
selection function applied to the EGS sample to balance 
the fraction of redshifts above and below z « 0.7. With 
the weight, Xii calculated in this way, we determine the 
differential galaxy stellar mass function: 

^{M^)dM^ = y -^dM, (6) 

j max 

4.1. Uncertainties and Cosmic Variance 

To estimate the uncertainty in our mass distribution 
we must account for several sources of random error and 
model their combined effect through Monte Carlo sim¬ 
ulations. Different error budgets are calculated for the 
primary spectroscopic and photo-z supplemented sam¬ 
ples via 1000 realizations of our data set in which we 
randomized the expected errors. In both samples we 
model uncertainties in Vmax arising from photometric er¬ 
rors, and simulate the error on the stellar mass estimates, 
which, as described in EH is typically 0.2 dex and is en¬ 
coded in the stellar mass probability distribution of each 
galaxy. 

For the primary spectroscopic sample, errors on the 
rest-frame {U — B) colors are estimated by noting the 
photometric errors for a given galaxy. We do not model 
the uncertainty in the [Oil] SFR because unaccounted 
systematic errors are likely to be greater than the ran¬ 
dom uncertainty of measurements of the [Oil] equivalent 
width. We stress again that these diagnostics are only 
used to separate the bimodal distribution into active and 
quiescent components. 

For the photo-z supplemented sample, we model red¬ 
shift uncertainties using Az/(1 + z) = 0.07 for the 
Rab < 24.1 ANNz subsample and Az/(1 + z) — 0.18 
for the Rab A 25.1 BPZ subsample. The BPZ uncer¬ 
tainty is slightly higher than measured in the comparison 
to the spectroscopic sample because we expect a poorer 
precision of BPZ on objects with Rab > 24.1. These 
redshift uncertainties affect the distribution of objects 
in our redshift bins as well as the Monte Carlo realiza¬ 
tions of stellar mass estimates in the photo-z sample. 
The final uncertainty at each data point in the stellar 
mass functions from both the spectroscopic and photo-z 
supplemented samples is estimated conservatively as the 
sum, in quadrature, of the Icr Monte Carlo errors and 
the Poisson errors. 

We now turn to cosmic variance. Because our sample 
is drawn from four independent fields, it is possible to 
estimate the effects of cosmic variance by comparing the 
results from different fields in each of our three redshift 
bins: 0.4 < z < 0.7, 0.75 < z < 1.0, and 1.0 < z < 1.4. 
For the spectroscopic sample with z > 0.75, we can com¬ 
pare the EGS to the sum of Fields 2, 3, and 4, yielding 
two subsamples with roughly equal numbers of galaxies. 
We compare the total and color-dependent stellar mass 
functions of these two subsamples and divide the median 
of the differences measured at all data points by \/2 to es¬ 


timate the cosmic variance. Unfortunately, galaxies with 
z < 0.75 come only from the EGS. The cosmic variance 
estimate here is derived by performing the same calcu¬ 
lation on three subsets of the EGS sample and dividing 
it by -s/S- Of course, cosmic variance on the scale of the 
EGS itself is not included in this estimate. For the three 
redshift intervals, this method provides average la sys¬ 
tematic cosmic variance uncertainties of 29%, 12%, and 
26%. 

In addition to these rough estimates, we have checked 
that the observed density distribution is not affected by 
cosmic variance. Based on the density-dependent mass 
functions from different fields, there is no evidence of a 
single structure or overdensity in one of the fields that 
would bias our results. We also note that the type- 
dependent mass functions are, to first order, affected by 
cosmic variance in the same way as the total mass func¬ 
tions. Thus, while absolute comparisons between differ¬ 
ent redshift intervals must account for cosmic variance 
errors, comparisons using the relative or fractional abun¬ 
dance of a given population are less susceptible to this 
uncertainty. For example, comparisons between different 
fields of the abundance of blue galaxies at z ~ 0.5 with 
masses near M* = IO^^'^Mq (see suggest that the ab¬ 
solute measurement is uncertain at the ~40% level, while 
the relative measurement made with respect to the to¬ 
tal abundance of galaxies is uncertain only at the ^15% 
level. For red galaxies in this same bin, the standard de¬ 
viation for the relative fraction is 10%, but increases to 
25% for the measurement of their absolute abundance. 

4.2. Completeness and Selection Effects 

We now turn to the important question of the redshift- 
dependent completeness of the stellar mass functions in 
our sample. Incompleteness resulting from the DEEP2 
selection technique is corrected for by applying weights 
to galaxies within the spectroscopic sample, as described 
in The R and Ag-band limits of this sample also 
introduce incompleteness effects, however, that cannot 
be corrected through weighting. Instead, we adopt two 
approaches to determine how the incompleteness from 
our magnitude limits affects the derived mass functions. 

First, we consider the likelihood of observing model 
galaxies of known mass and color based on the R and 
Ag-band magnitude limits of the primar y spectroscopic 
samp le. We follow previous work (e.g., iFontana et alJ 
l2f)nl and track the stellar mass of a template galaxy 
with a reasonable maximum M*/L ratio determined by 
models with solar metallicity, no dust, and a burst of 
star formation beginning at Zform = 5 that lasts for 0.5 
Gyr and is then truncated. We set the luminosity ac¬ 
cording to the R and Ag-band limiting magnitudes. The 
maximum stellar masses of these model galaxies provide 
estimates of our completeness limits. We consider two 
realizations. In the first, the model is placed at the low- 
z edge of each redshift interval where the intrinsically 
faintest galaxies would be detected (such galaxies have 
high Vmax weights). In the second, more conservative 
limit, the model is placed at the high-z edge of the in¬ 
terval, representing the faintest systems that could be 
detected across the entire redshift bin. The results indi¬ 
cate that the Ag-band limit ensures completeness above 
3 X IO^^Mq at all redshifts, while the i?-band limit pri¬ 
marily affects our high-z sample below IO^^Mq. We dis- 
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Fig. 3.— Apparent color-magnitude diagrams used to illustrate sample completeness. The plots show the distribution of {R — Ks) vs 
Ks for the primary spectroscopic redshift sample with Rab < 24.1 (solid, color circles) compared to the Rab < 25.1 sample, which has 
been supplemented with photometric redshifts (small black dots). The spectroscopic sample is colored according to location in the bimodal 
restframe {U — B) distribution as described in the text. The Rab = 24.1 and Rab = 25.1 magnitude limits are indicated by the solid 
diagonal lines. 


CUSS these limits in detail below. 

Our second approach is to measure the effects of in¬ 
completeness directly by comparing our primary spec¬ 
troscopic sample with Rab < 24.1 to the fainter sample 
with Rab < 25.1, supplemented by photometric red- 
shifts. This approach is particularly useful for investi¬ 
gating the way in which the i?-band limit introduces a 
bias against red galaxies, especially at z > 1. This bias 
could mimic the effect of downsizing by suppressing the 
fraction of red galaxies at z > 1. The behavior of model 
galaxies (described above) as well as the comparison to 
a fainter sample both yield consistent estimates for the 
sample completeness which we will show does not com¬ 
promise our results. 

Figure 0 compares the distribution in 2'.'0 diameter 
{R — Ks) vs total Ks color-magnitude space of the pri¬ 
mary spectroscopic sample with Rab < 24.1 (solid color 
circles) to that for the fainter Rab <25.1 sample (small 
black dots) supplemented with photometric redshifts. As 
expected, in the low redshift bin, the majority of the 
Rab < 25.1 sample is contained within the spectroscopic 


limit of Rab = 24.1 (although the photo-z sample in¬ 
cludes many more galaxies with Rab < 24.1 that were 
not selected for spectroscopy). At high redshift, however, 
the primary sample is clearly incomplete, with a substan¬ 
tial number of galaxies having Rab > 24.1. While the 
full range of (i? — Ks) colors is included in the primary 
sample, a color bias is introduced because the reddest 
galaxies are no longer detected at Kab ^ 20. 

To demonstrate how this color bias affects the mass 
completeness of the sample, the left-hand panels in Fig¬ 
ure 01 chart the stellar mass completeness of the photo-z 
supplemented sample limited to Rab < 24.1 as com¬ 
pared to the full sample with Rab < 25.1. With increas¬ 
ing redshift, the fraction of galaxies satisfying Rab < 
24.1 decreases. This results in a loss of detected galax¬ 
ies below a given mass which corresponds closely to the 
conservative limits based on model galaxies placed at the 
high-z edge of each redshift interval, as described above 
and indicated by the dotted vertical lines in Figure 0| 
However, because of the Vmax weighting, our analysis is 
actually more complete than the figure suggests. A more 
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Fig. 4. — Completeness of the mass distribution in the primary spectroscopic sample taking into account the R-band and Rs-band mag¬ 
nitude cuts. The left-hand panels show the completeness resulting from the Rab < 24.1 limit as derived from comparisons within the 
photo- 2 : supplemented sample. Vertical dotted lines indicate the mass of the faintest galaxy that could be detected at the hlgh-z end of 
each redshift interval, while vertical dashed lines correspond to galaxies at the low-z end. The right-hand panels, also drawn from the 
Rab ^ 25.1 sample, show the effect of the Rs-band limit on the mass completeness. Vertical dotted and dashed lines indicate model galaxy 
limits, as in the left panels. The more conservative dotted lines are adopted as the mass limits of the paper. 


appropriate completeness limit is indicated by the verti¬ 
cal dashed line, which corresponds to model galaxies at 
the low-z edge of each redshift interval where the faintest 
galaxies are still detected. 

Figure 01 indicates that for galaxies with z < 1, the 
i?-band limit has only a small effect on the spectroscopic 
sample. Galaxies with Rab < 24.1 constitute the vast 
majority of the total distribution in the first two redshift 
bins, with a minimum completeness of «80% in the mid¬ 
dle redshift bin. The i?-band limit has a greater effect 
in the high-z bin, where the Rab < 24.1 distribution 
accounts for just over «60% of the Rab < 25.1 sample 
over most of the mass range. We note that the i?-band 
limit affects mostly red galaxies which we will show dom¬ 
inate at higher masses. For this reason, the completeness 
function resulting the i?-band limit rises again at lower 
masses. In section we correct for the i?-band incom¬ 
pleteness and show that it does not affect our conclusions. 

The right-hand panels of Figure 01 illustrate the effect 
of the Ag-band limit on the mass completeness of the 
spectroscopic sample. Using our deepest Palomar obser¬ 


vations, we construct stellar mass distributions from the 
photo-z supplemented sample with Rab < 25.1. Again, 
for each redshift interval, we compare the mass distribu¬ 
tion of a subsample with a Ag-band limit equal to 0.5 
magnitudes fainter than the limit imposed on the spec¬ 
troscopic sample at that redshift, K^-^. As in the left- 
hand panels of Figure01 the resulting completeness func¬ 
tion yields an estimate of the Ag-band mass completeness 
limit that is consistent with the high-z and low-z model 
galaxy masses described above. For the analysis to fol¬ 
low, we conservatively adopt the high-z model limits indi¬ 
cated by the dotted lines, although we acknowledge that 
heavily obscured galaxies may still be missed above these 
limits. For each redshift bin, these limits correspond to 
10.1, 10.2, and 10.4 in units of log Mq. We note that 
for the middle and high redshift bins in which the R- 
band limit affects the sample, the Ag-band completeness 
limit occurs at lower masses, justifying our method for 
studying the R and Ag-band limits separately. 

5. RESULTS 
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5.1. The Total Stellar Mass Function 

The total galaxy stellar mass functions in three red- 
shift intervals are plotted in Figure 0 Results for both 
the spectroscopic, Rab < 24.1 sample (solid circles, dark 
shading) and Rab < 25.1, photo-z supplemented sam¬ 
ple (open diamonds, light shading) are presented. The 
width of the shaded curves corresponds to the hnal Icr er¬ 
rors using the Monte Carlo techniques discussed earlier. 
For the Rab < 24.1 sample, the lowest redshift interval 
draws only from the EGS, leading to a larger cosmic vari¬ 
ance uncertainty as indicated by the isolated error bar in 
the upper right portion of the plot. The Rab < 24.1 
mass functions in the two higher redshift bins, as well as 
the Rab < 25.1 sample in all bins, utilize data from all 
four DEEP2 survey helds. The close agreement between 
the spectroscopic and photo-z supplemented total mass 
functions indicates that the weighting applied to com¬ 
pensate for the DEEP2 spectroscopic selection function 
in the Rab <24.1 sample is successful at recovering a 
complete, magnitude-limited sample. The vertical dot¬ 
ted lines represent estimates of the mass completeness 
originating from the iC-band magnitude limit (see E3. 
The solid curve plotted for all three redshift intervals is 
the Schechter fit to the Rab < 25.1 mass function at 
z ^ 0.55, and can be used to gauge evolution within our 
sample. 

As discussed in S3 the mass functions of the primary 
sample presented in Figure |31 are affected by incomplete¬ 
ness because of the Rab = 24.1 spectroscopic limit. The 
degree of incompleteness is apparent in the comparison 
between the total mass functions of the spectroscopic 
and Rab < 25.1 samples. As expected from Figure 01 
the first two redshift bins are hardly affected by i?-band 
incompleteness. However, the Rab = 24.1 limit is clearly 
important in the high redshift bin where the comoving 
number density in the photo-z sample is larger by a fac¬ 
tor of «2 (0.3 dex) for M* < lOi°-®M 0 . 

Figure 01 also plots results from several previous stud¬ 
ies, all of which have been normalized to h = 0.7 and 
adjusted for comparison to the Chabrier IMF used here. 
The dotted curve is the Schechter fi t to the local stel - 
lar mass function as d etermined bv ICole et al.1 (1200111 . 
Results frorn the K20 (iFontana et 3,111200411. MUNICS 
llDrorv et aLll2004l . and COMB017 llBorch et alj|2006ll 
surve ys are also plotted. The K20 survey JCiinatti et all 
[ 200 l has good spectroscopic coverage (92%) but is shal¬ 
lower than the data presented here and covers only 52 
square arcminutes or 0.01 square degrees. The MUNICS 
survey llDrorv et al.ll200lll covers roughly one square de¬ 
gree, but only reaches K ^ 19.5 (Vega) and is composed 
of 90% photometric redshifts. The COMBO 17 results 
utilize carefully calibrated photo-z’s llWolf et all I2004)l 
and represent an area of 0.8 square degrees, but infrared 

S ' ' netry was not available. Instead, iBorch et all 
scale their stellar mass estimates to the magni¬ 
tude measured in a narrow filter (21 nm wide) at SIOA. 
This restricts the redshift range where masses can be es¬ 
timated to z < 1 . 

The stellar mass functions from previous work as well 
as the new results presented here are in good agree¬ 
ment in the first redshift bin in Figure 01 It is in¬ 
teresting to note that all of the results plotted hnd a 
higher abundance of massive galaxies (M* > 5 x IO^^Mq) 
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Fig. 5. — Total stellar mass functions in three redshift bins. 
Shading indicates the width of Icr error bars. The mass function 
for all galaxies in the primary spectroscopic sample is designated by 
solid circles with dark shading. That for the deeper, Rab ^ 25.1 
sample, supplemented by photo- 2 :’s, is indicated by op en diamonds 
and l ight gr ay shading. Results from the COMBOlTJBorch^^al] 
1200(31 . K20 IFontana et al .1200411 . and MUNICS iProrv et^l!n2004F 
surveys are also plotted. The dotted curve is t he Schechter fit to 
the local mass function from IHole et all I2001I1 . The solid curve 
is a Schechter fit to the Rab ^ 25.1 results in the first redshift 
bin. Vertical dotted lines show estimates for the mass complete¬ 
ness limit resulting from the i^s-band magnitude limits—in the 
highest redshift interval incompleteness resulting from the i?-band 
limit becomes important, as demonstrated by the difference be¬ 
tween the spectroscopic and photo- 2 : supplemented mass functions. 
The isolated error bar in the upper right-hand portion of each plot 
indicates the estimated average uncertainty due to cosmic variance. 


aud a lower abundance of less massive galaxies (M* < 
6 X 10^°M f7i ) com pared to the local mass functions of 
ICole et id] (1200111 . although a slight excess above the 
Schechter fit is a lso seen at the highest masses in lo¬ 
cal samples (e.g.. ICole et_alll^._ We do not plot the 
local mass function fronijBgll_g^^y because it is 

in goo d agreement with ICoDetafT l20f)]1l . IBorch et alJ 
(120001) argue that the higher number of massive galax¬ 
ies in their study may be caused by photo-z errors, but 
we find similar results in our spectroscopic sample. This 
discrepancy may reflect systematic differences in the way 
stellar masses are estimated, and highlights potential dif- 
hculties in comparing high-z work to local studies. 

Our mass function at z ~ 0.9 shows little evolution 
compared to the lower redshift interval. In fact, for 
M* > IO^^Mq, we find slightly higher number densities 
at z ~ 0.9, but this is likely due to cosmic variance. Our 
results are consistent with the COMBO 17 mass functions 
but are higher than the K20 and MUNICS mass func¬ 
tions, a difference that seems to increase in the highest 
redshift bin. Here we find only a slight decrease in the 
numbers of galaxies with M* > 3x IO^^Mq, an effect that 
is consistent with no evolution, given the uncertainties 
from cosmic variance. The photo-z supplemented mass 
function at z ~ 1.2 reveals a more significant decline. 
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TABLE 2 

SCHECHTER FiT PARAMETERS 


Redshift 

4>* 

Mpc“® dex“l) 

log hroM* 
(Mo) 

a 

Primary Spec- 2 :, Rab ^ 24.1 

0.4 < 2 < 0.7 
0.75 < 2 < 1.0 
1.0 < 2 < 1.4 

0.0027 ± 0.0004 
0.0031 ± 0.0010 
0.0012 ± 0.0008 

10.94 ± 0.08 
10.87 ±0.10 
10.97 ±0.25 

-0.81 ±0.2 

-0.59 ±0.4 
-0.51 ±0.5 

Photo -2 Supplemented, Rab ^ 25.1 

0.4 < 2 < 0.7 
0.75 < 2 < 1.0 
1.0 < 2 < 1.4 

0.0025 ± 0.0007 
0.0031 ± 0.0014 
0.0015 ± 0.0004 

10.93 ± 0.11 
10.87 ±0.10 
11.01 ± 0.07 

-0.95 ±0.1 
-0.51 ±0.4 
-0.93 ±0.1 


Note. — Quoted uncertainties were estimated based on field-to- 
field comparisons. 


however, with respect to our \ow-z bin in the number 
density of lower mass galaxies (with M* ^ 3 x IO^^Mq). 
This hints at the notion that the mass assembly history 
of all galaxies may proceeds in a top-down fashion, lend¬ 
ing support to rec ent claims of so-calle d “mass assembly 
downsizing” ('e.g.. iCimatti et ^l 200 ( 1 ^ . While intrigu¬ 
ing, we postpone a detailed analysis of this effect to a 
future paper, choosing instead to focus on downsizing in 
the context of star formation activity. 

In Table 13 we list the parameters of Schechter fits to 
the Rab < 24.1 and photo-z supplemented, Rab < 25.1 
mass functions plotted in Figure 0 By integrating the 
photo-z Schechter fits we obtain estimates for the global 
stellar mass density in units of log Mq /170 Mpc“^ of 
logp* = 8.31 ± 0.07 at z ^ 0.55, logp* = 8.30 ± 0.1 
at z ^ 0.9, and logp* = 8.15 ± 0.1 at z ^ 1.2, where 
the uncertainties come from cosmic variance estimates 
based on field-to-field comparisons. These results are 
in good agreement w i th previous measurements of Pt, 
llDickinson et alJl2003l: iRudnick et alJl2fiC)^ iDrorv et al.1 

1200 ^ 

5.2. The Mass Functions of Blue and Red Galaxies 

We now partition the total mass functions presented 
in FigureElinto active and quiescent populations accord¬ 
ing to the bimodality observed in the restframe {U — B) 
color. The results are illustrated in Figure 13 which for 
reference also plots the total mass functions from Figure 
13 As before, the vertical dotted lines indicate the onset 
of incompleteness resulting from the A^j-band limit, and 
the solid curve is a Schechter fit to the photo-z supple¬ 
mented total mass function in the first redshift interval. 

As discussed above, the R < 24.1 limit introduces in¬ 
completeness in the high redshift bin. To mitigate this ef¬ 
fect we derive a color-dependent completeness correction 
for the high-z bin based on the photo-z supplemented 
R < 25.1 sample. Inferring the restframe {U — B) color 
for this sample is difhcult because of photometric red¬ 
shift uncertainties. Instead, we adopt the simpler ap¬ 
proach of applying a color cut in observed {R — K), 
which, as shown in Figure |3 maps well onto the rest¬ 
frame {U — B) color for the high-z spectroscopic sample. 
We tune the {R — K) cut so that the resulting color- 
dependent mass functions of the Rab < 25.1 sample 
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Fig. 6.— Mass functions in three redshift bins partitioned by 
rest-frame (U — B) color as described in the text. Shading indicates 
the width of Icr error bars. The mass function for all galaxies in the 
primary spectroscopic sample is designated by solid circles. That 
for the deeper, Rab ^ 25.1 sample, supplemented by photo-z’s, is 
indicated by open diamonds. Vertical dotted lines show estimates 
for mass incompleteness resulting from the i^s-band magnitude 
limit. In the high-redshift bin (bottom panel), thin red and blue 
lines trace the expected increase in the red and blue populations 
if the spectroscopic limit of Rab = 24.1 were extended by one 
magnitude. The isolated error bar in the upper right-hand portion 
of each plot indicates the estimated systematic uncertainty due to 
cosmic variance. The solid curve in each panel shows the Schechter 
fit to the photo- 2 : supplemented mass function in the first redshift 
interval. 



Fig. 7.— Log fractional contribution of the red and blue popu¬ 
lations to the total stellar mass function at various redshifts. The 
relative contribution of red galaxies (left plot) increases with cos¬ 
mic time, while that for blue galaxies decreases (right plot). The 
expected high- 2 : fraction of red and blue types after correcting for 
R-h&nd completeness (as in Figure is shown by a thin white 
line for the red population and a thin blue line for the blue pop¬ 
ulation. The corrected fractions from the Rab ^ 25.1 sample are 
entirely consistent with those observed in the primary Rab ^ 24.1 
spectroscopic sample. Dashed lines in the right-hand panel show 
approximate fits to the declining fraction of blue galaxies. Mq is 
defined by where these lines cross the dotted horizontal lin e, wh ich 
marks a fractional abundance of 1/3 the total value, (see El. 
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match the spectroscopic {Rab < 24.1) mass functions 
above M* = where the spectroscopic sample 

is complete. This yields a value of {R — Ks) 2 " = 3.37, 
consistent with the (U — B) bimodality apparent in the 
color-magnitude diagram shown in Figure ^ 

Using this observed (R—K) color cut for the high-z bin 
only, we show the color-dependent mass functions for the 
Rab < 25.1 sample in Figure [Bias thin white and blue 
lines. While these curves suffer from their own uncertain¬ 
ties, such as photo-z errors and a less precise measure of 
color, they are useful for illustrating the nature of the 
i?-band incompleteness in the spectroscopic sample. It 
is important to note that while the deeper sample yields 
higher mass functions for both the red and blue popu¬ 
lations below « IO^^’^Mq, the relative contribution 
of each one to the total mass function is similar to what 
is observed in the spectroscopic sample. This is shown 
more clearly in Figure[71 It should also be noted that the 
Rab <25.1 sample itself is not complete, although the 
decreasing density of points near the Rab = 25.1 limit 
in Figure 0] suggests it is largely complete over the range 
of stellar masses probed. 

Figure El reveals several striking patterns. A clear 
trend is observed in which the abundance of massive blue 
galaxies declines substantially with cosmic time, with the 
remaining bulk of the actively star-forming population 
shifting to lower mass galaxies. As the abundance of the 
blue population declines, red galaxies, which dominate 
the highest masses at all redshifts, become increasingly 
prevalent at lower masses. The two populations seem 
therefore to exchange members so that the total number 
density of galaxies at a given stellar mass remains fixed. 
We also note the clear downward evolution of the cross¬ 
over or transitional mass, Mtr, where the mass functions 
of the two color populations intersect. Above Mtr, the 
mass function is composed of primarily red galaxies and 
below it blue galaxies dominate. We return to this be¬ 
havior in 

Figure [71 shows these results in a different way. Here, 
the log fractional contribution from the red and blue pop¬ 
ulations are plotted in the same panel so that the redshift 
evolution is clearer. The completeness-corrected color- 
dependent mass functions (with R < 25.1) are shown 
as the solid white and blue lines. Their overlap with 
the shaded curves from the high-z spectroscopic sample 
is remarkable and indicates that unlike absolute quanti¬ 
ties, relative comparisons between mass functions drawn 
from the spectroscopic sample are not strongly biased by 
the i?-band mass completeness limit. As noted previ¬ 
ously, plotting the relative fraction also removes the first 
order systematic uncertainty from cosmic variance, mak¬ 
ing comparisons across the redshift range more reliable. 

The downsizing evolution in Figure El is now more 
clearly apparent in Figure 0 The relative abundance 
of red galaxies with M* « 6 x increases by a 

factor of «3 from z ~ 1.2 to z ~ 0.55. At the same 
time, the abundance of blue, late-type galaxies, which 
are thought to have experienced recent star formation, 
declines significantly. 

5.3. Downsizing in Populations Defined by SFR and 
Morphology 

The patterns in Figure El and 0 are also apparent when 
the galaxy population is partitioned by other indicators 


of star formation. This is demonstrated in Figure0which 
shows the fractional contribution (in log units) of ac¬ 
tive and quiescent populations to the total mass func¬ 
tion. The “blue” and “red” samples dehned by restframe 
(U — B) color and shown in Figure El are reproduced here 
and indicated by solid circles. 

For the two redshift intervals with z > 0.75, we 
have plotted contributions from samples with high and 
low [OllJ-derived SFRs. We divide this sample at 0.2 
Mq yr“^, which is the median SFR of the star-forming 
population at 0.75 < z < 1.0. This imposes a more 
stringent criterion than the restframe {U — B) cut, which 
counts galaxies with only moderate or even recent star 
formation as “late-type.” Not surprisingly, the middle 
redshift bin contains fewer high-SFR galaxies compared 
to blue {U — B) systems and more low-SFR galaxies 
than red systems, although the mass-dependence ob¬ 
served with either criteria is qualitatively similar. 

In the high-z bin, the populations defined by color and 
[Oil] track each other more closely. Not only does this 
confirm that the mass-dependent evolution seen in Fig¬ 
ure El is reproduced when the sample is divided by the 
[Oil] SFR, but it also indicates that the average star for¬ 
mation rate is higher in this redshift bin. More of the 
“blue” population is now above the SFR cut as com¬ 
pared to the middle redshift bin. This evolution in the 
observed SFR will be discussed in detail in Noeske et al. 
(in preparation). 

It is helpful also to understand how these trends re¬ 
late to earlier work motivated by understanding the role 
of morphology in downsizing. Figure 0 also plots the 


contribution from morpholo gically-defined popu lations, 
drawing from the sample of iBimdv et in] il2n(^ which 
has been adjusted to the h = 0.7 cosmology used here (we 
note that the recent addition of HST/ACS imaging in the 
EGS provides an opportunity to ex tend this morpholog - 
ical comparison in the future). In iBundv et alJ l|2nn5jl . 
morphologies were determined visually using HST/ACS 
imag ing data from the GOODS fields llGiavalisco et alj 
O) and were divided into three broad classes: E/SO, 
spirals, and peculiars. The fractional contribution from 
the spiral and spiral-l-peculiar samples are plotted in Fig¬ 
ure 0 for comparison to the late-type populations de¬ 
scribed above. The E/SO fraction is compared to the 
early-types. I t should be noted that the smaller sam¬ 
ple size of the IBundv et alJ ll2005ll data leads to greater 
uncertainties and larger effects from cosmic variance. 

With these caveats, there is quite good agreement 
in the mass-dependent evolutionary trends between the 
morphological and color/SER selected samples. In de¬ 
tail, the fraction of ellipticals is systematically lower 
than the red/low-SFR populations while the fraction 
of spirals-bpeculiars is systematically higher than the 
blue/high-SFR galaxies. This suggests that the process 
that quenches star formation and transforms late-types 
into early-types operates on a longer timescale for mor¬ 
phology than it does for color or SFR. We return to this 
point in m We also note that spiral galaxies do not 
always exhibit star formation and can be reddened by 
dust while som e ellipticals have experienced recent star 
formation ('e.g.. lTreu et al.l200,5ll that could lead to bluer 
colors. 
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Fig. 8. — Fractional contribution in log units to the total mass function from quiescent (left-hand plots) and active populations (right- 
hand plots) defined by rest-frame {U — B) color, [Oil] star formation rate, and morphology (see text for det ails). Star formation rates 
based on [Oil] are only available for galaxies with 2 : > 0.75. The morphological mass functions are those of IBundv et al ] J2003). The 
i^s-band completeness limits are indicated as in Figure]^ 


5.4. Quantifying Downsizing: the Quenching Mass 
Threshold, Mq 

Several authors have identified a characteristic tran¬ 
sition mass, Mtr, which divides the galaxy stellar mass 
function into a high-mass regime in which early-type, 
quiescent galaxies are dominant and a low-mass regime 
in which late-type, act i ve galaxies aue dominant (e.g., 

Kju^mann et al.l2003tilBaldrv et al.l2004fBundv et al'l 

200 , 5tl . Using our various criteria, the downward evo¬ 
lution of this transitional mass with time is clearly 
demonstrated in the upper panel of Figure For the 
color based mass functions, good agreement is found 
in comparison to the COMBO 17 results discussed in 
lBorch_et_alJ12006 ). The morphological sample is taken 
from IBiindv et al' II2005I1 . where we have grouped spi¬ 
rals and peculiars into one star-forming population and 
compared its evolution to E/SOs. In the high-z bin, Mtr 


for the morphological sample occurs beyond the probed 
stellar mass range and so we have extrapolated to higher 
masses to estimate its value (this uncertainty is reflected 
in the horizontal error bar at this data point). 

The color-defined Mtr shows a redshift dependence of 
Mtr oc {1 + 2 )“^, similar to that for the morphological 
sample. Stronger evolution is seen for the the [OII]- 
defined samples as expected if evolution is more rapid for 
the most active sources. As discussed in the mass 
scale of morphological evolution is approximately 3 times 
larger (ftiO.5 dex) than that defined by color or [OB]. We 
also note from FigureElthat Mtr does not change appre¬ 
ciably when the i?-band mass incompleteness is corrected 
in the high-z bin. 

While the evolution in Mtr is illustrative of downsiz¬ 
ing, since its definition in terms of physical processes is 
completely arbitrary (equality in the relative mass contri- 
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Fig. 9.— Redshift evolution of the characteristic transitional 
mass, Mtr (top plot), as well as the quenching mass, Mq (bottom 
plot). The plotted redshift values indicate the approximate average 
redshift in each bin. In both plots, the observed behavior is shown 
for the population partitioned according to morphology (light, open 
squares), restframe (U — B) color (solid circles), and [Oil] SFR 
(asterisks). 

butions of two populations), its significance is not clear. 
We prefer to seek a quantity that clearly describes the 
physical evolution taking place. Accordingly, we intro¬ 
duce and define a quenching mass limit, Mq, as that 
mass above which star formation is suppressed in galax¬ 
ies. This threshold is a direct byproduct of the mecha¬ 
nism that drives downsizing. We define Mq by noting 
the mass at which a line fit to the declining fraction of 
star-forming galaxies (Figure [7|l drops below 1/3. We 
apply this simple definition to the other measures of star¬ 
forming populations and plot the resulting values of Mq 
in the bottom panel of Figure El and list them in Tablejol 
The relative behavior of differently classified populations 
is similar to the top panel, but the physical mass scale 
associated with Mq is a factor of 2-3 higher than that 
of Mtr- We find an approximate redshift dependence of 
Mq (X {z + 1)'^'®, similar to the dependence of Mtr- 

Mq is a useful quantity because it reveals the masses at 
which quenching operates effectively as a function of red¬ 
shift. What is more, blue and star-forming galaxies at the 
quenching mass also tend to be among the reddest and 
least star-forming as compared to the full late-type popu¬ 
lation as a whole at each redshift, suggesting such galax¬ 
ies are in the process of becoming early-types. As we 
discuss below, the quantitative evolution of Mq strongly 
constrains what mechanism or mechanisms quench star 
formation and also provides a convenient metric for test¬ 
ing galaxy formation models. 

5.5. The Environmental Dependence of Downsizing 

We have so far considered the mass-dependent evolu¬ 
tion of late- and early-type populations integrated over 
the full range of environments probed by the DEEP2 
Redshift Survey. We now divide the sample by environ¬ 
mental density to investigate how this evolution depends 
on the local environment. 

In Figure El we plot galaxy stellar mass functions for 
the samples shown in Figure El split into low-density en¬ 
vironments on the left-hand side and high-density envi- 


TABLE 3 

The Quenching Mass Threshold, Mq 


Redshift 

Blue [U - B) 

High SFR 

Spirals+Peculiars 

0.4 < 2 ; < 0.7 
0.75 < ^ < 1.0 
1.0 < ^ < 1.4 

10.73 ±0.08 
10.98 ±0.09 
11.23 ±0.13 

10.80 ±0.10 
11.25 ± 0.18 

11.20 ± 0.16 

11.65 ±0.29 
12.00 ±0.55 


Note. — Mq is given in units of log(Mfe^Q/MQ). Morphology 
data comes from|Bund^_e^_^ | 2003 ). 


ronments on the right-hand side. The density discrimi¬ 
nant is simply the median density (indicated by an over¬ 
density of zero in Figure EJ, so all galaxies are included. 
The evolution of the total mass functions in the two den¬ 
sity regimes follows similar patterns as for the combined 
sample in Figure El Furthermore, the relative contribu¬ 
tion of active and quiescent galaxies, as quantified by 
Mtr and Mq, exhibit the same kind of downsizing signal 
seen in Figure El Crucially, there is little difference (less 
than 25%) in the value of Mtr between the below- and 
above-average density samples. 

The absence of a strong environmental trend in the 
downsizing pattern of massive galaxies is a surprising re¬ 
sult that warrants further scrutiny. Although no strong 
de pendence was seen in the Fundamental Plane analysis 
of lTreii et alJ ll2005ll . their density estimates were much 
coarser and based on photometric redshi f ts. Lo cal work 
utilizing the SDSS by iKaiiffmann et al.1 ll2004ll demon¬ 
strated differences in galaxy properties as a function of 
environment, but covered a much larger dynamic range 
in density than is accessible with our sample. We empha¬ 
size that the results of this work focus primarily on field 
galaxies, and, with the typical volumes probed (~10® 
Mpc^), only rich groups are included at the highest den¬ 
sities. DEEP2 is too small to probe cluster scales. 

It is important to note that individual density measure¬ 
ments, even in our spectroscopic analysis, can exhibit 
uncertainties of a factor of 3, so there is some overlap 
expected between the above- and below-average density 
regimes in Figure El However, because the total mass 
functions exhibit different abundances—with a factor of 
^2 more high-mass galaxies in the above-average density 
regime—the above- and below-average density subsam¬ 
ples successfully probe different populations. The higher 
abundance of massive galaxies at increa sing densities 
agrees well with trends in local studies (e.g.,|Balog hjetal| 
l200lL[20?)^lBaldrv et al.ir2004 IKaiiffmann et al.ll2004ll . 

Still, the comparison above is inevitably dominated by 
galaxies near the median density (Figure EJ, and so we 
adopt a second, more stringent approach that divides 
the distribution into three, more extreme density bins. 
In this way, we can sample the full dynamic range of our 
survey. As described in H3.3I we define a second set of 
density thresholds at ±0.5 dex from the median density 
and construct mass functions for the extreme high- and 
low-density regimes. We plot the results in Figure ITTl 
which again shows the fractional contribution of the red 
and blue populations to the total mass function. This 
time, the samples are also divided by density, with the 
high density points connected by solid lines and the low 
density points connected by dotted lines. 
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Fig. 10. — Stellar mass functions in the three redshift bins at above- and below-average density partitioned by rest-frame {U — B) color. 
The vertical dotted lines show the i^s-band magnitude completeness limits. The estimated cosmic variance is designated by the isolated 
error bar. Significant incompleteness from the i2-band limit is expected in the high redshift bin, but it is not possible to apply corrections 
as in Figure 13 because the local density can only be measured in the spectroscopic sample. The solid line shows the Schechter fit to the 
photo -2 supplemented total mass function in our lowest redshift bin. 


Figure m illustrates a more substantial environmen¬ 
tal effect in the mass-dependent evolution of the sample 
that in each redshift i nterval is suggestive of l ocal envi¬ 
ronmental trends te.g. ^Kauffmann et al.ll200'3l . Our ob¬ 
servations allow us to investigate how these trends evolve 
with time. We find that the rise of the quiescent pop¬ 
ulation and the evolution of Mtr appears to be acceler¬ 
ated in regions of high environmental density. This effect 
does not depend on the particular choice of the density 
threshold, although the differences between the two en¬ 
vironments grow as more of the sample near the median 
density is excluded from the analysis. The difference in 
Mtr between these two environments is roughly a factor 
of 2-3. We caution that interpreting the results in the 
high-z bin is difficult because the effects of completeness 
and weighting are most important here. However, the 
fact that the Mtr values in the two density regimes are 
more similar in the high-z bin (an offset of ~0.1 dex) 
might suggest that the structural development or physi¬ 
cal processes that lead to the density dependence at lower 
redshifts are not fully in place before z ^ 1. 


The environmental dependence observed in Figure ITTI 
does not occur because of the high-density regime alone, 
as one might expect given the potential presence of dense 
structures in the sample. The environmental effect is less 
strong but still apparent in comparisons between the high 
and middle density regimes as well as between the middle 
and low density regimes. 

In summary, our various measures of downsizing, Mtr 
and Mq, depend strongly on redshift (they evolve by 
factors of ~3-5 across the full redshift range), but less 
so on environmental density for most field galaxies (less 
than 25%). In support of previous inferenc es on the role 
of environment based on local studies (e.g..lBalogh et alJ 
12004 iKaiiffmann et al1l2004 iHla.nton et al.ll2004l . this 
direct measurement of density-dependent evolution is 
particularly striking and serves to emphasize that galaxy 
mass, not environmental location, is the primary parame¬ 
ter governing the suppression of star formation and hence 
producing the signature of downsizing. 


6 . DISCUSSION 
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Fig. 11.— Fractional contribution in log units to the total mass 
function. In each of three redshift bins solid lines connect points 
representing the extreme high-end of the density distribution, while 
dotted lines indicate the extreme low-end. The Ks-band complete¬ 
ness limits are indicated as vertical dotted lines. 


6 . 1 . The Rise of Massive Quiescent Galaxies 

Figures El and 0 demonstrate a clear feature of the 
downsizing signal observed since z ^ 1 , namely the in¬ 
crease in the number density of massive quiescent galax¬ 
ies. Although our results are consistent with previous 
studies which found a rise in the red g alaxy abundanc e 
of a factor of ~ 2 - 6 , depending on mass llBell et aDiOOl, 
the present work represents a significant step forward not 
only in its statistical significance and precision by virtue 
of access to the large spectroscopic and infrared data set, 
but also in clearly defining the mass-dependent trends. 

In discussing our results we begin by considering 
the processes that might explain the present-day pop¬ 
ulation of early-type galaxies. In order to reconcile 
the significant ages of their stellar populati ons implied 
by pr ecise Fundamental Plane studies (e.g., iTreu et ^ 
l2nnl E^iTder We] et a,l 11200,^ wit h hierar c hical modek 
of structure formatio n, Bell et af] ll2005aD . iFaber et al.1 
l|200,^ . Ivan DokkumI l|200,^ and others have introduced 
the interesting possibility of “dry mergers”—assembly 
preferentially progressing via mergers of quiescent sub¬ 
units. 

While dry merg ers clearly occur ijTran et alJ l200fil: 
lyan DokkumI 1200,^ . our results suggest that they can¬ 
not be a substantial ingredient in the assembly history 
of massiye quiescent galaxies. As shown in Figure El 
the obseryed increase in the number of quiescent sys¬ 
tems is almost perfectly mirrored by a decline in star¬ 
forming galaxies such that the total mass function ex¬ 
hibits little eyolution over 0.4 < z < 1.4 compared to 
that of the red and blue populations. A simple transfor¬ 
mation of one into the other is sufficient to high preci¬ 
sion. For example, simply interchanging the numbers of 
red and blue galaxies in the high-z bin in the mass range 
1011.2 ^ logM*/M 0 < 10 ^^ ®, leads to a prediction for 


the number density in the middle-z bin that is accurate 
to within ^25%, well within the cosmic variance uncer¬ 
tainty. In addition, galaxies defined as blue or late-type 
with masses near Mq tend to be among the reddest and 
least star-forming, suggesting they are likely candidates 
for becoming early-type systems. 

It is conceivable that the dry merger rate is mass- 
dependent and conspires to move galaxies along the mass 
function in a way that leaves its shape preserved. This 
would imply the presence of massive galaxies at low red- 
shift that are not seen in our sample. For example, 
using the approx imate dry merger rate estimate from 
iBell et al.l ll2005aD of 1.3 mergers every 6.3 Gyr, ^25% 
of red galaxies would have to experience a major dry 
merger (in reality, this rate, calculated over 0 < z < 0.7, 
might be expected to be higher at z ~ 1 ) between 
the high-z and middle-z bins. If we apply this rate 
of assembly to the red population in our high-z bin at 
log(M*/MQ) = 11.4, we find that the total mass func¬ 
tion in the middle-z bin at log(M*/M 0 ) = 11.7 should 
be higher by approximately 60% compared to what is ob¬ 
served. While this one data point represents only a «2cr 
result, similar arguments apply across the mass function 
and between the middle-z and low-z bins. More detailed 
modeling of the effect of merger rates on the mass func¬ 
tion will be presented in a further paper (Bundy et ah, 
in preparation). 

Most difficult for the dry merger hypothesis as a key 
feature of galaxy formation is the weak environmental 
dependence we observe in the downsizing signal. Given 
a mechanism has to be found to preferentially bring qui¬ 
escent sub-units together, one would expect a strong den¬ 
sity dependence in the dry merger rate. By contrast, in 
our large sample it is clear that the majority of quiescent 
galaxies were assembled in a manner more sensitive to 
mass than environment. 

Finally, our results show that downsizing is not only a 
feature of the star formation histories of massive quies¬ 
cent galaxies (e.g., ITreu et af.llW)^ | van der Wei et alJ 
uni, but is also apparent in the way their abun¬ 
dance increases. This observation is of particular inter¬ 
est for galaxy formation models based on the hierarchi¬ 
cal ACDM framework. By incorporating AGN feedback 
effects into the semi-analyti c models from the Munich 
group, iDe Lucia et alJ (1200511 predict star formation his¬ 
tories for massive ellipticals that follow the downsizing 
trend in which more massive galaxies host older stel- 
lar populations . How ever, as shown in Figure 5 from 
IDe Lucia et al ] (l 2 f)f)l . these models still predict a hierar¬ 
chical mass assembly history for e llipticals. The p r oblem 
in comparing the predictions of IDe Lucia erall l|2005ll 
to our results (Fi gure [71 for example) is that the simu¬ 
lated ellipticals in IPe Lucia et al.l (1200511 are defined by 
their appearance at z = 0 , not the redshift at which 
we observe them. While the observations presented here 
suggest a top-down pattern in the growing abundance of 
early-types, it is not yet clear how significant a problem 
this may be for current semi-analytic models. 

6 . 2 . The Origin of Downsizing 

The results of this paper reveal important clues as to 
the nature of downsizing and, via a clear measurement 
of the trends, will assist in constraining and ruling out 
several of the popular explanations. A detailed compari- 
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Fig. 12. — Fractional contribution in log units to the total mass 
function, divided into stellar mass bins and plotted as a function of 
redshift. The plotted redshift values indicate the approximate av¬ 
erage redshift in each bin. The vertical dashed line in the top panel 
indicates the onset of incompleteness from the /?-band limit which 
occurs above the redshift indicated. The if^s-band completeness 
limit does not affect the results plotted in this figure. 


son with such models is beyond the scope of the current 
paper but we discuss some of the key issues here. 

First we wish to dismiss a possible suggestion that our 
discovery of the quenching mass threshold is somehow an 
artifact of our selection process. For example, it might 
be argued that our result could arise from a uniform de¬ 
crease in the incidence of star formation at all masses 
combined with a survey selection effect in which rare, 
massive objects are seen only at higher redshifts because 
of the larger volumes probed. Our results exclude this 
possibility. First, our sampled volumes are relatively sim¬ 
ilar (only a factor of ~4 difference between low-z and 
high-z), and we demonstrate that the fraction of star¬ 
forming galaxies depends not only on redshift but mass 
as well. Figure IT^ illustrates this point. Here we divide 
the sample into smaller redshift intervals and follow the 
evolution of the red and blue populations in the three 
largest stellar mass bins, charting the fractional contri¬ 
bution of the two populations to the combined mass func¬ 
tion. The highest mass bin contains the largest fraction 
of quiescent galaxies at almost every redshift, and the 
transformation of the active population into passively- 
evolving systems occurs first in the high mass bin and 
later at lower masses. The rate of change in the inci¬ 
dence of star formation is clearly mass-dependent. 

Turning now to physical explanations, we have shown 
in H5.4l that the quenching mass threshold, parametrized 
by Mq, provides a very useful description of how the 
fraction of star-forming galaxies evolves. The question 
then is what mechanism is responsible for this quench¬ 
ing? Can it adequately reproduce the quantitative trends 
observed, for example the environmental dependence? 
Merging may provide a starting point for answering this 


question and explaining the transformation of late-types 
into early-types. Merging between disk systems has long 
been thought to be an important mechanism by which 
ellipt i cals form le.g. lToonirelll977t iBarnes & Hernaiii^ 
11 9911 iSnringel et al.ll2nn5bll and the similar behavior of 
the morphological and color-defined values of Mq in 
Figure M suggests that the same process governing the 
growth of ellipticals may also broadly explain the rise of 
quiescent galaxies in general. However, as argued above, 
significant merging is likely to affect the shape of the total 
mass function which does not appear to evolve strongly 
in our sample. Merging may be occurring at masses be¬ 
low our completeness limit, but the observed evolution in 
the relative mix of early- and late-type galaxies suggests a 
process that quenches star formation first with morpho¬ 
logical transformation following later. This would lead 
to higher values for Mq in the morphological samples, as 
observed. 

Merger-triggered quenching has further difficulties. 
Fundamentally, the hierarchical merging of dark matter 
halos is expected to proceed from low mass to high mass, 
not the other way around. One solution to this difficulty 
would be to appeal to the fact that merging and assem- 
bly rates are accelera ted in regions of high density (e.g., 
iDe Lucia et al ] |200l which also host the most massive 
systems. Over a range of environments, downsizing could 
arise naturally from the fact that massive galaxies live 
in these accelerated environments. However, we find no 
significant density dependence in the bulk of our sam¬ 
ple (Figure [TOIl . with downsizing appearing in all envi¬ 
ronments. This suggests that density-dependent merger 
rates are not the answer and that an internal feedback 
process on galactic scales is largely responsible for driv¬ 
ing the downsizing pattern. 

Many groups have recently suggested that internal 
AGN feedback may be the missing ingredient. Trig¬ 
gered perhaps by merging, energy deposited by the AGN 
effectively quenches further star formation, eventually 
transforming late -type galaxies into early-types (e.g., 
iSilk fc Heeill199^ . There are at least two current im¬ 
plementations of this process. In the explosive model, 
energy from an intense “quasar phase” (usually triggered 
by merging) heats and expels gas in the h alo that would 
otherwise be availabl e for sta r formation llGranat O-eli-^ 
12004 iSoringel et alJ l2005at iScannaoieco et alJ 1200,^ . 
While inspired by hydrodynamic simulations, it is not 
clear how the downsizing mass scale could arise from 
this model or what prevents future cooling of the ex¬ 
pelled gas. In the “radio feedback” model, low luminos¬ 
ity AGN energy couples to hot gas in halos with masses 
greater than Mohock in which the cooling time is longe r 
than the dynamical time (see lDekel fc Birnboimll^004ll . 
This results in a long-term suppression of star forma¬ 
tion in massiv e halos and, on even larger scales, cool¬ 
ing flows (e .g..lCrqton et al.ll2005t IDe Lucia et al.ll2005t 
iBower et all 1200,4 iGattaneo et all l200(ill . These mo dels 
produce older stellar populations in more massive galax¬ 
ies but still predict that low mass early-types (as defined 
at z = 0) assemble first (see Figure IT)). While a possi¬ 
ble mass scale is provided in Mohock that may help drive 
the downsizing signal, the physical mechanism that cou¬ 
ples the radio feedback to the hot gas is not well under¬ 
stand. Beyond simulations and formation models incor¬ 
porating AGN, additional evidence for the role of AGN 
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feedback comes from phenomenological studies that re¬ 
late observations of black hole mass, the quasar luminos¬ 
ity function, an d the properties of galaxies in a consistent 


framework (e. 

Ef.. IMerloni et al.l200^ 

llHopkins et al.l2005l: 

IMerloni et abl 

l2nn6HHonkins et abl 

2 ()n^. 


Regardless of the physical explanation for downsizing 
in the star formation activity of galaxies (and several 
may be necessary), it is clear that precise quantitative 
measures of the evolving mass distribution and its depen¬ 
dence on the basic parameters explored here will provide 
the ultimate test of these theories. 

6.3. Reconciling Downsizing with the Hierarchical 
Structure Formation 

In this section we step back from specific models of 
galaxy formation and explore in a more general sense 
the kinds of processes that may reconcile downsizing 
with hierarchical structure formation. Over the past 
decade there has been strong confirmation from many 
independent observations that the large-scale structure 
of the universe matches the basic predictions of the cold 
dark matter (CDM) model in which dark matter “halos” 
harbor galaxies and grow through constant, hierarchical 
merging with other halos. The masses of all halos in¬ 
crease monotonically with time, but at every epoch, it 
is the largest halos that are growing most actively and 
that are also the last to have assembled. Thus at hrst ex¬ 
amination, downsizing, in which star formation proceeds 
from high-mass to low-mass systems, seem at odds with 
the CDM picture. 

Several processes may be contributing to reversing the 
bottom-up trend in CDM structure formation to pro¬ 
duce what appears to be a top-down pattern in galaxy 
formation. The first is a gradual effect resulting from 
the accelerating expansion of the universe which causes 
halo growth rates to slow once the universe reaches a 
scale factor (1-1-2;)“^ > D,m- The second is the physics 
of gas cooling, which selects a characteristic mass scale 
for galaxy formation llRees fc Ostrikedll977t lSilklll97^ 
IWhite fc ReeslIfOT^ . Gas cannot cool rapidly, and by 
implication stars cannot form efficiently, until structure 
formation produces halos with virial temperatures in ex¬ 
cess of 10"^ K . This sets the epoch for the onset of galaxy 
formation at 2; ~ 15-20. Once gas temperatures in¬ 
crease to 10®-10^ K, cooling once again becomes inef¬ 
ficient, turning off star formation in the most massive 
halos. This then marks the end of galaxy formation, as 
more and more mass builds up in group and cluster halos 
above this cooling limit. 

More detailed numerical or semi-analytic models of 
galaxy formation show that the cooling delay alone is 
insufficient to reduce star formation to observed levels, 
particularly in massive halos, and that other forms of 
feedback are required. While the details remain contro¬ 
versial (e.g., iBensou et 3.1.1120?)^ . the net effect of this 
feedback is to place an upper limit on the range of halo 
masses over which active star formation can take place. 
This limit helps explain why star formation in galaxies is 
rarer at the present-day than it was at z ^1-2. However, 
as discussed in Hti.2l it is less obvious how to explain the 
observed decline in the mass scale of star-forming galax¬ 
ies. 

To help gain insight into this question, we can attempt 
to relate various galaxies in our sample to dark mat¬ 


ter halos. Models of the halo occupation distribution 
function predict that galaxies with the range of stellar 
masses sampled here {\og{M/M q) ~ 10-12, correspond¬ 
ing to logLfej ~ 10-11.2) should res ide in dark mat¬ 
ter halos of mass lO^^-lO^^M fr^ fe.g.. lYang et al.ll200,3l: 
iCoorav fc Milosavlievi(i2005a|) . and furthermore that 75- 
80% of these galaxies will be “central,” that is the dom- 
inant galaxies within their hal o, rather than satellites 
ifCoorav fc Milosavlievicll2005bll . These models suggest 
that the objects in the three mass bins in Figure^l corre¬ 
spond approximately to central galaxies in galaxy, group, 
and cluster halos. 

In the top panel of Figure El we show the comov¬ 
ing number density of halos of mass \og{M/M q) = 

12.5.13.5 and 14.5 as a function of redshift (three lines, 
from top to bottom). The numbers are roughly consis¬ 
tent with the comoving number densities of galaxies in 
our three mass bins, although the most massive stellar 
objects are more abundant than Mq halos and may 
therefore reside in slightly less massive systems. The bot¬ 
tom panel shows the mean ages of halos in the three mass 
bins as a function of observed redshift. The mean age 
here is defined as the time elapsed since half of the halos 
in that mass range had hrst built up 90% or 50% of their 
current mass in a single progenitor (solid and dashed 
curves respect i vely), calculated using equation 2.26 from 
iLacev fc C^ (Il993|). 

By comparing Figure cni to Figure d we can gain in¬ 
sight into the relationship between the assembly of dark 
matter halos and the process that drives the observed 
mass-dependent decline of star-forming galaxies. The 
fraction of the most massive blue galaxies declines from 
nearly 100% to ^15% by z « 1 in Figured For the 
least massive galaxies, this same decline lasts over the 
full redshift range of our sample, implying a timescale 

2.5 times longer. Coupled with Figure El~which in¬ 
dicates that the least massive halos in our sample are 
about twice the age of the most massive halos—this sug¬ 
gests a process that, for the low-mass systems, is at least 
5 times slower than the growth of the dark matter halos 
in which they reside. Since global dynamical timescales 
should be independent of halo mass at a given redshift, 
this suggests that the quenching mechanism is strongly 
mass-dependent with the potential for different physical 
processes acting in different mass ranges. 

7. CONCLUSIONS 

Using a large sample that combines spectroscopy from 
the DEEP2 Galaxy Redshift Survey with panoramic 
near-IR imaging from Palomar Observatory, we have in¬ 
vestigated the mass-dependent evolution of held galaxies 
over 0.4 < z < 1.4. We have constructed stellar mass 
functions for active and quiescent populations, dehned in 
several ways, and divided into different samples accord¬ 
ing to accurate measures of the environmental density 
determined from the extensive spectroscopic data. We 
summarize our conclusions below: 

• The mass functions of active and quiescent galax¬ 
ies integrated over all environments conclusively 
demonstrate a downsizing signal. We quantify this 
by charting the evolution in a “quenching mass,” 
Mq, which describes the mass scale above which 
feedback processes suppress star formation in mas- 
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Fig. 13.— (Upper panel) The abundance of halos like ly t o host 
central galaxies in the three mass ranges plotted in Figure^] versus 
redshift. The curves are labeled with log(M/M 0 ). (Lower panel) 
Mean age of these systems as a function of their observed redshift. 
The mean age is defined as the time elapsed since half of the sys¬ 
tems had built up a fraction / of the mass they have at Zobs • Solid 
curves show the age for / = 0.9 and dashed curves that for / = 0.5. 


sive galaxies. We find that Mq cx (1 + z)^'^ with 
a factor of wS decrease across the redshift range 
probed. 

• The growth in the abundance of quiescent or 
“early-type” galaxies occurs first at the highest 
masses and then proceeds to lower mass systems. 
The relative abundance of early-types with M* ~ 
10 ^°Mq has increased by a factor of ~3 from z ~ 
1.2 to z ^ 0.55, whereas the total mass function ex¬ 
hibits little evolution in shape and normalization 
(less than 0.2-0.3 dex). This implies that early- 
type systems result largely via the transformation 
of active star-forming galaxies, indicating that “dry 
mergers” are not a major feature of their assembly 
history. 

• Alternative ways of dividing active and quies¬ 
cent galaxies, including the use of [Oil] equiva¬ 
lent widths and HST-derived morphologies, show 
qualitatively similar mass-dependent evolution and 
quenching. Interestingly, we observe that morpho¬ 
logical evolution appears to take place on longer 
timescales than changes in the apparent star for¬ 
mation rate which operate at lower mass scales at 
each redshift. 

• For the majority of galaxies in our sample, down¬ 
sizing shows little dependence on environment. An 
environmental signal is apparent when the ex¬ 
tremes of the density distribution are compared. 
In this case, downsizing in high-density regimes 
appears moderately accelerated compared to low- 
density ones, with values of Mtr lower by a factor 
of ~2. 


• We discuss several possibilities for the origin of 
downsizing based on our results. We clearly rule 
out a scenario in which the incidence of star forma¬ 
tion decreases uniformly for galaxies at all masses. 
The weak density dependence also argues against 
explanations that rely on the accelerated assembly 
of structure in dense environments, favoring inter¬ 
nal mechanisms instead. 


• Through comparisons to the expected behavior of 
dark matter halos, we argue that the dynamical 
timescale resulting from the growth of structure is 
at least 5 times longer in galaxies hosted by halos 
with log M/Mq « 12.5 (log M*/Mq « 10.8) com¬ 
pared to logM/MQ = 14.5 (logM*/M0 11.6). 
Because global dynamical scales are also indepen¬ 
dent of halo mass at a given redshift, this suggests 
that the quenching mechanism is strongly mass- 
dependent with the potential for different physical 
processes acting in different mass ranges. 


There are two obvious avenues for further studies of 
downsizing. In a forthcoming paper (Bundy et ah, in 
preparation) we discuss the constraints on merging and 
the growth of galaxies determined by our observations of 
the total mass function. This will help dissect the role of 
merging in driving downsizing. In the near future, it will 
also be possible to chart the incidence of AGN among 
the galaxy population and compare it to the incidence 
of star formation to probe the link between quenching 
and AGN. The significant Chandra follow up observa¬ 
tions currently underway in the EGS will make that field 
particularly exciting for such work. Other efforts from 
the DEEP2 collabora tion have provided new insight on 
environmental trends ijGooner et al.ll2?ir)^ and will focus 
on precise measures of the evolving star formation rate 
(Noeske et ah, in preparation). 
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